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ABSTRACT 

Wang, Quan. Ph.D., Purdue University, May 2012. Charge Multiplicity Asymmetry 
Correlation Study Searching for Local Parity Violation at RHIC for STAR Collabo- 
ration . Major Professor: Fuqiang Wang. 

The strong force is one of the four fundamental interactions in particle physics de- 
scribing the interaction between partons (quarks and gluons) which make up hadrons. 
The theory of the strong force is called quantum chromodynamics (QCD), which is a 
quantum field theory of the color charged partons. The force between color charges 
does not diminish while they are separated. This property causes the color charges 
to be confined in to hadrons, in ordinary matter. Quark-Gluon Plasma (QGP) is 
one phase of the QCD matter at extremely high temperature and/or pressure, where 
the partons are asymptotically free. Experimentally, QGP might be created in ultra 
relativistic heavy ion collisions [IH3]. 

It has been suggested that in such deconfined QCD matter, the metastable do- 
mains with non-zero topological charge Q w will generate charge separation along the 
system angular momentum direction caused by chiral magnetic effect (CME). The 
charge separation direction is random as the sign of Q w is random from domain to 
domain. The event-by-event charge separation along the system angular momentum 
direction violates the parity and time-reversal symmetries locally (LPV) [51410] . In 
this analysis, we measure the CME/LPV in heavy ion collisions with charge multi- 
plicity asymmetry correlations. 

We separate a heavy ion collision event into up and down, or left and right 
hemispheres according to the reconstructed event-plane and the plane perpendicu- 
lar to the event-plane. We then calculate the multiplicity asymmetries of the posi- 
tive and negative charges by taking the multiplicity difference between up and down 
hemispheres (A± >UD ), as well as left and right hemispheres (A± >LR ), and divide by 
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the total multiplicities. Since the event-plane does not distinguish between up and 
down nor left and right, the average of the asymmetries are consistent with zero, 
(A±,ud) = {A± t Ln) = 0. However, the correlations between the asymmetries are 
non-zero due to the physical correlations between particles. 

We study the variances ((A^ D ) and (A 2 LR )) and covariances ((A + A_) UD and 
(A + AJ)lr) of the charge multiplicity asymmetries. The asymmetries are calculated 
using the multiplicity from one side of the TPC tracks with respect to the event-plane 
reconstructed from the other side of the TPC tracks in order to avoid self-correlation. 
We also apply single particle detector efficiency correction on asymmetry calculation 
and event-plane reconstruction. The variance results are alway positive because they 
are the square of real numbers, which is the effect of statistical fluctuations. We 
subtract the statistical fluctuation and the detector non- uniformity effects to obtain 
the dynamical variances 5(A 2 ID ) and 5(A 2 LR ). 

We show the dynamical variances and covariances of Au+Au 200 GeV collisions 
and d+Au 200 GeV collisions. Data show the dynamical variances are positive at 
peripheral collisions consistent with d+Au data. This suggests that same-sign particle 
pairs are emitted preferentially in the same direction. Both variances in UD and LR 
drop in mid-central and central collisions and turn to negative, which suggests that 
the same-sign pairs are more likely to be emitted symmetrically, more back-to-back 
in other words, regardless of the directions. The covariances are largely positive for 
both UD and LR directions, which suggests the opposite-sign particles are strongly 
correlated, and emitted with small angle correlation. 

The CME/LPV expects charge separation across the event-plane (UD direction), 
which gives additional correlation to the same-sign particle pairs in out-of-plane (UD) 
direction than in-plane (LR) direction, i.e. a wider distribution of the asymmetries 
in UD direction. We should expect A(A 2 ) = 5(A 2 JD ) — S(A 2 LR ) > 0. One also expects 
that the positive and negative charges are anti-correlated in UD direction, so that the 
covariances A(A + A_) = (A + A_)ud — (A + A_) lr are negative. We show the UD — 
LR correlations of the dynamical variances and covariances. Both the variance and 
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covariance differences are positive for all centralities except the most peripheral bins. 
The variance UD — LR correlation is positive, which is consistent with CME/LPV 
expectation. However we also know same-sign pairs are preferentially back-to-back 
from mid-central to central collisions. The covariance UD — LR correlation is also 
positive, which is not consistent with the naive expectation of CME/LPV. 

We study the transverse momentum (jpx) dependence of the correlations. The 
CME/LPV expects the charge separation is mostly a low-p^ effect. However data 
show the correlations increase with in the mid-central collisions. 

The dynamical variances and covariance as a function of event-by-event anisotropy 
v^ s for low-pr (pt < 2 GeV/c) and high-pr (p? > 2 GeV/c) particles are analysed. 
The UD — LR variance and covariance show opposite trend of low-pr vf" 3 , but with 
very weak dependence of the high-p^ v^ s . We use four different data and cuts to 
verify the result: sub-events with r/ > and rj < 0, sub-events with large pseudo- 
rapidity gap 7] > 0.5 and 77 < —0.5, events with the first order ZDC-SMD event-plane 
and top 2% most central data. They all show similar v^ s dependence. 

There might be charge independent common background sitting between the same- 
sign and opposite-sign correlations. So we define charge separation as the difference 
of same-sign and opposite-sign correlations, A = A(A 2 ) — A(A + A^) , to cancel the 
background, and show it as a function of the wedge size, azimuthal region of the 
analysed particles. The charge separation vanishes with the decrease of the wedge 
size, which suggests the charge separation effect is within the vicinity of the reaction- 
plane. 

We also show the charge separation as a function of event-by-event low-p^ v^ s . 
The CME/LPV effect does not expect event anisotropy dependence. However, we 
see the charge separation is strongly and linearly depending on \ow-px v^" 5 for four 
different cases. The linear dependence intercept at zero or slightly positive, when the 
sub-event is isotropic in low-pr particle azimuth angle, i.e. v^ s = 0. 

Because the average v^ s is positive due to elliptic flow, the charge separation is 
positive if we integrate over all events. The linear v^ s dependence of charge separation 
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is more likely an intrinsic charge dependent bulk correlation of the medium. A more 
precise measurement of the CME/LPV effect may lie on the events with a more 
isotropic shape. We show the charge separation of events with \v^ s \ < 0.04 as a 
function of centrality, which is consistent with zero. Thus we give the upper limit of 
the CME/LPV effect based on the linear fit of the charge separation A as a function 
of v° bs , that A = 4.7 x 10~ 5 with 98% CL. 
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1. INTRODUCTION 

I. 1 Strong Interaction and Quantum Chromodynamics 

There are four fundamental interactive forces in particle physics, which describe 
the way elementary particles interact with each other. They are electromagnetism, 
strong interaction, weak interaction and gravitation. Modern physics attempts to 
explain all physically observed phenomena by the theories of these fundamental inter- 
actions. Except for gravitation, theories of electromagnetism, the strong interaction 
and the weak interaction are well established in the standard model. In the concept 
model, matter consists of elementary particles, which are spin one-half fermions and 
interact with one another according to their properties (charges). They interact by 
exchanging spin one gauge bosons, also called force carriers. A summary of the fun- 
damental interactions with their theories and properties is shown in table 11.11 Figure 

II. II shows the three generations of elementary leptons and quarks, as well as the gauge 
bosons. 

Specifically in this thesis, the focus is on the study of the strong interaction. It 
is a short range interaction comparing to other three interactions that binds protons 
and neutrons together to form the nucleus of atoms (in the range of 1-3 fm), and also 
binds quarks and anti-quarks to form hadron particles (in the range of less than 1 
fm). The strong interaction is carried out by gluons exchanging "color charge", an 
analogous to the electronic charge in electromagnetism between quarks, anti-quarks 
and gluons. Different from the force carrier photon (7) in electromagnetism, gluon (g) 
can interact between themselves. Unlike electromagnetism' s electric charges (positive 
and negative), there are three types of color charges, resulting in different behaviors 
of strong interaction. The behavior of the color charges and the interactions of quark- 
gluon are detailed in the theory of quantum chromodynamics (QCD), the quantum 
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Table 1.1 

Fundamental interactions. 



Interaction 


Theory 


Force Carriers 
(gauge boson) 


Relative 
Strength 


Strong 


Quantum 
ChromoDynamics (QCD) 


gluon (g) 


10 38 


Electro- 
magnetic 


Quantum ElectroDynamics 
(QED) 


photons (7) 


10 36 


Weak 


Electroweak Theory 


W (W*) and 
Z (Z°) bosons 


10 25 


Gravitation 


General Relativity (GR) 


gravitons (hypo- 
thetical) 


1 



field theory of the nuclear interaction. As part of the standard model, mathematically, 
the theory is a non-Abelian gauge theory based on a local symmetry group SU(3). 

There are two unique properties of the QCD theory: color confinement and asymp- 
totic freedom. Quarks and gluons are the only elementary particles carrying color 
charges. The strong force between color charges, unlike all other forces, doesn't di- 
minish with increasing the distance of the color charges. It takes an infinite amount 
of energy to separate two quarks. Thus, before the quarks can be separated, the 
energy is large enough to create quark and anti-quark pairs to combine with the 
original quarks. Experimentally, isolated quarks have never been observed, i.e. free 
color charges. Any ordinary matter which can be observed are color neutral. This 
phenomenon is called color confinement. 

Asymptotic freedom is a property of the SU(3) gauge theory. At high energy, or 
equivalently at very short distance (<^C 1 fm), the interaction between quarks becomes 
weak, while at low energy or equivalently at large distance (> 1 fm), the interaction 
becomes strong. This phenomenon prevents the baryons and mesons from unbinding. 
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Figure 1.1. The three generations of quarks and leptons, and the gauge 
bosons. Figure is taken from wikipedia. 



1.2 Quark-Gluon Plasma and Chiral Symmetry 

Quark-Gluon Plasma (QGP) is a phase of QCD matter at extreme conditions, 
such as very high temperature (T) and/or high baryon chemical potential hb [TT|IT2"]. 
Figure 11.21 shows an illustrated phase diagram of quark matter. QGP might be 
produced in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) [IH1] , 
which is similar to the environment of the universe a few milliseconds after the Big 
Bang. The hot dense matter created in the heavy ion collisions at very high energy 
is considered to be in thermal equilibrium. Further studies show that extra degrees 
of freedom could be released at the quark level. Quarks have been relatively freed 
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from the confined nuclei. This phenomenon is called color deconffnement, one of two 
fundamental properties of QGP. 



Figure 1.2. The illustration phase diagram of quark matter. 

The other fundamental property, also considered cLS cL SI gnature of the QGP, is 
chiral symmetry restoration. Chirality in physics is a phenomenon that an object 
is not identical to its mirror image, i.e. its mirror image cannot be mapped by 
only rotations and translations. In the high energy limit, chirality can be treated as 
helicity (handedness). The helicity is defined as the sign of projecting the particle spin 
onto its direction of motion. The chirality is positive (right-handed) if the direction 
of the particle's spin is aligned in the same direction as its motion. It is negative 
(left-handed) if the spin is opposite to its motion. 

For massless particles, the helicity cannot be reversed by a Lorentz boost because 
no observer can travel faster than light. Therefore the massless particles have their 
helicity fixed for all reference frames, such as photon (7) and gluon (g). Their helicity 
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are invariant under Lorentz transformation. On the other hand, massive particles 
have slower speed than light. Then one can always boost a reference frame, so that 
the momentum reverses the direction. Massive particles thus may change their helicity 
signs after the Lorentz boost. 

In QCD theory, the Lagrangian can be written as 

/ / 

where / and g denote flavor index and the QCD coupling constant, F^ u denotes the 
spin-1 gluonic field strength tensor. is the vector potential of the color field, ipj 
and t a are the quark fields and the generators of the color SU (3) group. The mass 
term mj^fipf explicitly breaks the chiral symmetry of the QCD Lagrangian. 

Quark masses come from two sources. One is the "naked" quark mass, also called 
current mass, which is considered to originate from the Higgs mechanism in standard 
model. The other source is from the gluon field induced by a valence quark (quark 
which determines the hadron's quantum number), where the quark is surrounded like 
a cloud by sea quarks called covering. The two terms together give rise to the effective 
quark mass called the constituent mass. 

For light quarks, i.e. up and down, the constituent mass is much larger than the 
current mass, while for heavy quarks, i.e. charm, bottom and top, the constituent 
mass is nearly the same as the current mass. Figure 11.11 shows the current masses 
of the quarks. To show the constituent quark mass, we use proton as an example. 
The proton is a composite of three valence quarks ct uud" with the total current mass 
approximately 10 MeV. However, the mass of proton is 938 MeV, which is much 
larger than the total current masses of the valence quarks. The difference comes from 
the gluon field, the binding energy of quantum chromodynamics, while the gluons are 
massless. 

In the QGP phase as shown in the phase diagram, the quarks can be released 
from the confined matters and move relatively "freely". They are not really free, 
but relatively free. Thus, in such quark matter, quarks loose their covering. The 
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light quarks will have their mass greatly reduced to nearly zero, such that the mass 
term in the QCD Lagrangian vanishes, which results in the chiral symmetry becoming 
restored in the quark matter. In such chiral limits (m u = 0, = 0), for light quarks, 
all left-handed quarks remain left-handed, and all right-handed quarks remain right- 
handed. Each chiral state has a chiral symmetry partner with the opposite parity 
and equal mass. 

1.3 Chiral Magnetic Effect and Local Parity Violation 

Based on the well defined gauge theory, many remarkable properties of QCD 
matter have been discovered. One of the properties is that the field configurations 
can be characterized by a topological invariant, the topological charge Q w [10]. It is 
defined as 



gluonic field tensor and its dual. 

In QCD matter with the chiral limits (massless quarks rrif = 0) satisfied, chiral 
symmetry can be restored in the initial state with Nlh — N^h, where Nm and 
Nrh denote number of left-handed and right-handed quarks. There are metastable 
domains with certain topological charge Q w forming in the vicinity of the deconfined 
QCD matter [6], which lead to parity (V) and charge-conjugation and parity (CV) 
violation, if Q w is non-zero. The result is to convert right-handed (left-handed) quarks 
to left-handed (right-handed) quarks, with Nlh ~ Nrh — 2A r jQ w in the final state, 
where Nf is the number of flavors [TU] . 

Figure 11.31 shows an illustration of the Chiral Magnetic Effect (CME) inside the 
quark matter with the presence of a large and uniform magnetic field B. Within the 
quark matter, all the quarks are deconfined, and chiral symmetry is restored with 
chiral limit (m u = = 0). The red arrows denote the momentum direction, and the 
blue arrows denote the spin direction of the quarks. Due to the large magnetic field, 




(1.2) 



where g is the QCD coupling constant, and F" and F^ u 



= \e^ v p(T F^ denote the 



B 







Figure 1.3. Illustration of the chiral magnetic effect in a very large homo- 
geneous magnetic field. The red arrows denote the momentum direction, 
and the blue arrows denote the spin of quarks. See text for the details of 
the effect. Figure is taken from reference [TO] . 



quarks will eventually occupy the lowest Landau level after equilibrium, with their 
magnetic moments align in the same direction as the B field. Thus, u quarks with 
positive charges have their spin in the same direction of the B field, and d quarks with 
negative charges in the opposite direction of the B field as shown in part (1). Then 
the gauge field with non-zero topological charge Q w interacts within the metastable 
domain, as shown in part (2), and breaks the chiral symmetry by, for example with 
negative Q w , converting left-handed u and d quarks into right-handed quarks. The 
result is to flip the momentum direction of the left-handed quarks to the opposite 
direction. In the end as shown in part (3), all the u quarks are right-handed and 
moving upwards carrying positive charges, and all the d quarks are also right-handed 
but moving downwards carrying negative charges. This charge separation effect is 
then called the Chiral Magnetic Effect (CME), and could possibly be measured in 
experiment if indeed true. 



8 



parity violation 




(a) Parity violation (b) Charge separation 



Figure 1.4. Panel (a): Illustration of parity violation with charge sepa- 
ration in the system angular momentum direction. Panel (b): Illustration 
of charge separation along the system angular momentum direction (L) 
and magnetic field direction (B) in center of mass frame. The z direction 
is defined as the projectile (A) and target (B) nuclei momentum direction. 
The two heavy ion nuclei undergo a mid-central collision with the reaction 
area shown in orange color. See text for details. 



The schematic view of the charge separation effect in a mid-central heavy ion 



collision is illustrated in figure 1.4(b) In the center of mass frame, the beam direction 
(z direction in the figure) and the direction connecting the centers of two colliding 
nuclei (x direction) define the "Reaction-Plane" of the collision. The overlapping 
reaction area, has an elliptical shape and contains the hot dense QCD matter which 
could have chiral symmetry restored. The spectators, the wounded nuclei A and B 
in the figure, carry positive charges and create a magnetic field B when passing the 
reaction area. The direction of the magnetic field is the same as the QCD system's 
angular momentum L direction. If CME is indeed true and the effect can survive 
through the hot dense medium evolution to the detectors, one should observe charge 
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separation along the magnetic field, i.e. system angular momentum direction, which 
is perpendicular to the reaction plane. 



As shown in figure 1.4(a), charge separation gives the system an electric dipole 
moment (EDM), with its direction pointing from the negative charge to the posi- 
tive charge. When the system has its angular momentum aligned (anti-aligned) with 
the direction of the EDM, applying a parity operation to the system will change its 
parity state to anti-alignment (alignment). Thus, charge separation in the angular 
momentum direction is a phenomenon of parity violation. It is known that the par- 
ity symmetry is violated in some weak interactions, while it is well preserved in all 
other three interactions including strong interaction. If we could measure the charge 
separation in the system angular momentum direction in QCD matter, it may indi- 
cate that the parity symmetry could be violated in strong interaction. However, the 
topological charge Q w of the metastable domains are random within the QCD mat- 
ter. Thus, the direction of charge separation is also random. The charge separation 
effect cancels out over repeating experiments, in other words, it only happens on the 
event-by-event basis. So, the effect is only local, namely local parity violation (LPV). 
Globally, the parity symmetry is still conserved over a large number of events. 

The estimates of the charge separation signal are proportional to the topological 
charge Q w , and diluted by the event multiplicity. It was first calculated in reference [7] 
that the asymmetry of u quarks, for example, can be estimated as 

when assuming V and CP are preserved in the hadronic process. Then the asymmetry 
should translate to the hadron multiplicity asymmetry of charged pions 

A,^f + , (1.4) 

where iV + is taken as positively charged pion multiplicity in one unit of rapidity. 
This is because soft particles are usually correlated in one unit of rapidity range, 
which is also the extent of the parity violation domain in the rapidity space [?]. At 
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STAR experiment, the reference multiplicity NR e fMuit is the total number of charged 
particles recorded by the main TPC with pseudo-rapidity range of —0.5 < r] < 0.5, 
see section 12.2.41 At RHIC energy of 200 GeV per nucleon pair in Au+Au central 
collision, N^ e fMuit is typically around 300, and drops to 150 in mid-central collisions, 
and 30-50 in peripheral collisions [13J. Although the gold nucleon carries positive 
charges causing positive charged particles to produced slightly more frequently than 
the negatived charge particles (about 2%), we still have A + = A_ because they are 
normalized by the total number. Putting those numbers together, we have the charge 
multiplicity asymmetries at the order of ~ 10~ 2 in mid-central collisions [7HS], and 
~ 10~ 4 to ~ 10~ 3 for the asymmetry correlations. 

Apparently we do not take into account some of factors in the above estimation 
which may vary the final result, such as the magnetic field strength and duration time. 
More accurate estimations can be found in reference [10J, where theoretical calcula- 
tions have included magnetic field strength and fluctuations. They all give similar 
estimated results. However, the in-medium interaction and final state interactions 
before freeze out may also play an important role. Because the strong interaction 
conserves parity symmetry, the in-medium interactions could not contribute to the 
parity odd signal, but they will destroy the charge separation signal which are gen- 
erated in early stage of the collisions. The effect is to smear out the correlations 
between charges. The estimated effect is to reduce the signal for at least one order of 
magnitude [H], with the estimate around ~ 10~ 5 to ~ 10~ 4 for the asymmetry corre- 
lations [10]. A recent estimate shows that the CME/LPV induced charge asymmetry 
correlation is less than 10~ 6 in [T5], and 10 -4 after multiplied by the N part . Some even 
claim it is possible that the radial flow can even push the opposite-sign pairs into the 
same direction [16]. Thus, the sign of the opposite-sign correlation may possibly be 
even positive. 
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1.4 Anisotropic Flow and Three-Particle Correlator 

To study the properties of the QCD matter, physicists study heavy ion collisions 
at ultra relativistic conditions by accelerating heavy nuclei such as gold to the speed 
close to the light, and colliding them to create the new state of hot dense matter. A 
massive amount of particles are created during the collision, which mimics the early 
time of the Big Bang of the universe. 



In non-central heavy ion collisions shown in figure 1.4(b) , the initial spacial anisotropy 
will cause a pressure gradient in the azimuthal angle. The pressure in in-plane direc- 
tion is larger than that in out-of-plane, which translates to larger transverse momen- 
tum (jo?) and more particles are emitted in-plane than out-of-plane in the final state 
hadrons. The spacial and momentum space anisotropy of the event is driven by the 
initial pressure gradient, and affected by the medium strong interactions, which can 
be used as a probe of the initial collision geometry and the medium properties. The 
anisotropy is characterized by the Fourier expansion of the event azimuthal angle 0: 

dN 



d0 



oc 1 + ^ 2v n cos [n ((f) - i(>rp)] , (1.5) 



n=l 



where tjjRp denotes the true reaction plane as shown in figure 1.4(b) The Fourier 
coefficient v n stands for the n-th harmonic of the event azimuthal anisotropy. If we 
apply orthogonal condition of cosn(</> — ippp) to the above equation, we can get the 
n-th harmonic coefficient v n as: 

v n = (cos[n((p - 4>rp)]), (1.6) 

where (...) denotes an average over all the particles of each event. Due to the reflection 
symmetry, the sine terms vanish. Specifically, we refer the first order harmonic V\ 
as directed flow, and the second order harmonic Vi as elliptic flow, respectively. The 
directed flow v\ is an odd function of the rapidity due to momentum conservation, 
which is measured very small at RHIC y/s^N = 200 GeV Au+Au collisions with 
typically \vi\ < 0.005 for |^| < 1 [17] . The elliptic flow Vi is an even function of the 
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rapidity, and is measured sizable positive at 200 GeV Au+Au mid-central collisions 
about 6 percent for particle within 0.15 < px < 2.0 GeV/c [IS] . 

Both the first and second order anisotropy are directly related to the initial con- 
dition of the collision. Thus, the directed flow and elliptic flow provide us the exper- 
imental tools to determine the reaction-plane direction. The methods are introduced 
in sections 13.5.11 and 13.5.21 

Equation 11.51 is parity even because cosine is an even function to the mirror reflec- 
tion. In order to study the P-violation across the reaction-plane, we have to account 
for the parity odd terms, sine. The modified Fourier expansion can be written as 

d N °° 

— - oc 1 + 2 (v n cos [n ((f) - tpRp)} + a n sin [n ((f) - V>ap)]) • (1.7) 
d(p ^ 

n=l 

The coefficient a n stands for the "P-violation terms across the reaction-plane. As 
we introduced in previous section, a n is due to the local parity violation with the 
topological charge Q w . The signs of a n vary with the fluctuation of Q w . If we average 
a large amount of events, the averages of a n vanish because the topological charge is 
random. Thus, the direct measurement of a n is not possible. However, the charge 
separation effect will not vanish, and can be measured through correlation methods. 

STAR has published measurements of the first order charge dependent coefficient 
ai correlations. A charge dependent three-particle correlator [TI) |[T9~ll2T)] is introduced 
as 

(cos(0 a + (f) P - 24> c ))/v 2<c » (cos(0 Q + (pp - 2i) RP )) (1.8) 

= (cos A(f> a cos A^^) — (sin A(f> a sin A(f> b ) (1.9) 
= [{vi, a vifi) + B in ] - [(ai ;Q ,a 2 ,/?) + B out ], (1.10) 

where a, (3 and c are particle charge labels, and A0 = (f> — ipRp refers to the particle 
azimuthal angle relative to the reaction-plane. Assuming firstly, directed flow term 
(vi, a Vi,p) vanishes because it is an odd function of the rapidity and its fluctuation 
is small. Secondly, the average background from in-plane Bi n and out-of-plane B out 
cancels out, assuming the reaction-plane dependent background [B in — B out ] is small. 
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Lastly, only the first order a± dominates. Under such assumptions, the three-particle 
correlators are reported as the first evidence of the CME/LPV. We will review the 
assumptions and compare our observables to the three-particle correlations in section 
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2. EXPERIMENT 

In this chapter, we introduce the experiment of relativistic heavy ion collision. We 
introduce the facility and detectors used for data taking. We also introduce the 
kinematic variables measured for this analysis. 

2.1 Relativistic Heavy Ion Collider 

The Relativistic Heavy Ion Collider (RHIC) is one of the high energy heavy- 
ion colliders located at Brookhaven National Laboratory in Upton, New York on 
Long Island. By accelerating and colliding heavy ion and polarized proton beams, 
physicists study the matter created at extremely high temperature and density, which 
is the QCD matter with strongly interacting partons (quarks and gluons). 

Protons and heavy ion nuclei are accelerated in two independent pipes to nearly 
the speed of light and may collide at four intersecting points where the pipes cross. 
So far, several particle species have been accelerated for collisions at different en- 
ergies, including proton-on-proton (p+p), deuterium-on-gold (d+Au), copper-on- 
copper (Cu+Cu), gold-on-gold (Au+Au), copper-on-gold (Cu+Au) and uranium- 
on-uranium (U+U). For heavy ion collisions, the center of mass energy can reach 200 
GeV per nucleon pair. For p+p collisions, it achieved 500 GeV in 2009. 

As shown in figure |2~TT| the RHIC accelerator ring is 3,834 m long in circumference, 
and there are four experiments on RHIC collision points. They are STAR (6 o'clock), 
PHENIX (8 o'clock), PHOBOS (10 o'clock) and BRAHMS (1 o'clock). PHOBOS 
and BRAHMS have completed their commissioning and been shut down after 2005 
and 2006, while STAR and PHENIX are still running since 2000. This thesis is based 
on the data taken by STAR experiment of Au+Au and d+Au collisions at ^/saT/v = 
200 GeV. 



1(3 




Figure 2.1. The schematic plot of the Relativistic Heavy Ion Collider 
(RHIC) facility located at Brookhaven National Laboratory (BNL). Dia- 
gram taken from |21j . 



2.2 STAR Experiment 

The Solenoidal Tracker at RHIC (STAR) detector is located at the 6 o'clock 
interaction region of the RHIC accelerator ring. The main physics goal is to study 
the formation, evolution and characteristics of the strongly coupled Quark Gluon 
Plasma (sQGP) [IH1], a state of QCD matter which is believed to be formed at very 
high temperature and/or high energy density. It is designed primarily for charged 
hadron production measurements with high precision tracking and momentum over 
a large solid angle. 
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FTPC I 

Figure 2.2. Schematic STAR detector layout with selected sub-systems. 
Diagram taken from [21 j . 

STAR is a massive detector weighing 1,200 tons and as big as a building. It has 
a large relatively uniform acceptance in azimuthal angle, and also a large polar angle 
coverage in mid-rapidity. It consists of several sub-systems as shown in figure 12.21 
System upgrades are constantly going on since the first build of the detector. Some of 
the sub-detectors are removed and replaced for better measurement or more tracking 
capability. 

For convenience, we set up the STAR coordinate system with its origin located 
at the STAR detector geometry center. The z direction points to the west, which is 
parallel to the beam pipe direction. And the x direction points to the south, and the 
y direction points to up. The whole detector is surrounded by the main magnetic coil 
which generates a field in the z direction with the maximum of \B Z \ = 0.5 T |22j . 
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2.2.1 Time Projection Chamber 

The main tracking device is the Time Projection Chamber (TPC), which records 
the tracks of particles and provides the kinematic information of each track J22JIM]- 
It is located within the magnetic coil, with a 4.2 m cylinder length, 0.5 m inner radii 
and 2 m outer radii. With the full magnetic field (\B Z \ = 0.5 T) turned on, the 
TPC can identify a broad transverse momentum range of charged particles from 0.15 
GeV/c to 30 GeV/c depending on particles. 




Figure 2.3. The schematic cutaway view of STAR TPC detector. Figure 
is taken from [24] . 

A diagram of the TPC is shown in figure [2731 The TPC is a particle detector which 
consists of a gas-filled cylindrical chamber with multi-wire proportional chambers on 
the endcaps at each side of the cylinder. There is a high voltage central membrane 
disc which divides the cylindrical chamber into two halves. The central membrane, 
together with the "Outer Field Cage" and "Inner Field Cage" as shown in figure I2.3[ 
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provide a nearly uniform electric field along the z direction parallel to the beam pipe 
and magnetic field, pointing from the endcaps to the center. 

When a charged particle is generated in the collision, it traverses the TPC volume, 
ionizes gas atoms every few tenths of a millimeter along its trajectory, and leaves a 
cluster of electrons behind. The electron clusters then drift within the electric field 
toward the sectors of the endcaps. Each electron cluster will be accelerated by the 
electric field around the anode and cause a localized cascade of ionization, which is 
collected on the high voltage wire and results in an electric current. The position of 
the hit point is then recorded in the r and dimension, and the amplitude of the 
current is proportional to the energy of the detected particle. The endcaps are made 
by 24 identical sectors, 12 on each side. Each sector covers about 7r/6 in azimuthal 
angle, and totally full coverage of 2n. The z position is obtained by measuring the 
drift time from the collision to the time when the electron cluster is recorded at the 
endcaps because the drift velocity can be precisely measured beforehand. One can 
then reconstruct the tracks by fitting the 3-dimensional hit points collected by the 
TPC. Track in the TPC usually has a helix shape, and the maximum hit points at the 
endcaps can be as many as 45 hits. The transverse momentum pt and the charge sign 
of a particle can be calculated from the curvature of the trajectory with the applied 
magnetic field. Figure |2~41 shows the side view (xz plane), and front view (xy plane) 
of an Au+Au collision event. 

In this thesis, the particle information is taken mostly from the TPC. The most 
frequently used kinematic variables and 77 are shown in figure 12.41 The azimuthal 



angle (figure 2.4(a) ) is the angle of a particle transverse momentum p? relative to 
x axis in STAR coordinates. Pseudo-rapidity rj describes the angle relative to the 
beam direction, i.e. z direction. It is defined as 



77 = — In 



tan 
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(b) Side view 



Figure 2.4. An example event of central Au+Au collisions in transverse 
view (a) and side view (b). The definition of azimuthal angle <fi and 
pseudo-rapidity t] are shown in panel (a) and (b) respectively. 



where 6 is the angle of the particle momentum p relative to the beam axis, and pl 
denotes the longitudinal component of the particle momentum p. In the relativistic 
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limit when particle's speed is close to the speed of light, or the particle mass is small 
compared to its total energy, the rest mass could be ignored to good approximation. 
Then pseudo-rapidity is numerically close to rapidity which is defined as 



where E is the energy of the particle. The rapidity definition is used in theoretical 
calculations, while it requires two parameters E and p of a particle to be measured at 
the same time, which is inconvenient experimentally. However, the pseudo-rapidity 
requires only one measurement, p, and it is close to rapidity for light particles (i.e. pi- 
ons, electrons). The STAR detector has a large and uniform pseudo-rapidity coverage 
of -1 < 77 < 1. 

When an event is recorded, the tracks reconstructed from the TPC hits are called 
global tracks. By extrapolating the tracks back to the center of the detector, the 
original collision vertex can then be fitted from all the trajectories very precisely. 
The position of the collision vertex is noted as v in STAR coordinates. We can 
then do the tracking reconstruction again with all the TPC hits plus the additional 
collision vertex. The tracks reconstructed with collision vertex are called primary 
tracks. The advantage of primary tracks is that particles created from the collision 
vertex will get a better resolution since the collision vertex has a very good resolution. 
However, the tracks from secondary decay will have worse resolution because the 
collision vertex might be off from the track helix. Thus, we introduce the Distance of 
Closest Approach (DCA). It is defined as the closest distance from the collision vertex 
to a track helix. By cutting on DCA, one can manually select the tracks preferentially 
from the collision vertex or the secondary decay vertex. 

2.2.2 Zero Degree Calorimeters 

The Zero Degree Calorimeter (ZDC) detectors are located at ±18 m away from 
the geometry center of the STAR detector as shown in red in figure 12.21 [25] . They 
measure the energy deposition of neutrons, which are associated with the spectator 




(2.3) 
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matter, in the three tungsten plates on each side of the ZDC The ZDC has been used 
for beam monitoring and event triggering. 

There was an upgrade in 2003 of the ZDC by adding Shower Maximum Detectors 
(SMD), which gives the ZDC-SMD the capability of recording the shower profile 
of neutron clusters in the transverse plane (xy plane) of the collisions. The SMD 
information contains 7-slate (vertical) by 8-slate (horizontal) readouts from the ADC 
photomultiplier tubes which connect to the SMD scintillators. The raw readouts are 
corrected for the background by subtracting the pedestal. The readouts also have 
to be adjusted for the distortion of the electronics by applying the gain correction, 
which is from the cosmic ray calibration. Finally, the vertical and horizontal signals 
can present a well defined spectator position in the xy plane, which can be used for 
the direct event-plane reconstruction. 

2.2.3 Event Triggering 

The collision rate at RHIC 2004 is about 10M Hz, which means there are 10M 
collision events per second happening in the STAR detector. However, the STAR 
TPC is a slow detector because the electrons need time to drift to the endcaps to be 
recorded, and the DAQ system is also limited by the bandwidth. At 2004, the DAQ 
operated at rates about 100 Hz of Au+Au 200 GeV collisions. 

In order to reduce the recorded collision events, trigger detectors are needed to 
select 100 events out of 10M events per second, based on our interests [26]. The 
fast detectors are used as trigger detectors. They are the Zero Degree Calorimeters 
(ZDC), the Central Trigger Barrel (CTB) and the Beam-Beam Counters (BBC). The 
ZDC is introduced in previous section. The CTB detector is located between the 
outer cage of the TPC and the ToF detector, which measures the charge multiplicity 
in the same pseudo-rapidity range of |^| < 1 as the TPC with full azimuthal coverage. 
The BBC detector is located outside of the pole tip magnets as shown in figure I2T21 It 
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measures the multiplicity in forward region and provides vertex location information 
of the collision. 

The Au+Au minimum bias collisions, the least biased data sample, are triggered 
by both ZDCs (east- and west-side) and the CTB. The event is cut on the ZDC 
coincidence rate above certain value from the east and west ZDC, and the primary 
vertex from the ZDC signals. And there is also a cut on the CTB multiplicity to reject 
the non-hadronic events. However, some of the low multiplicity hadronic events are 
also rejected, thus there are bias on the low multiplicity events at the most peripheral 
collisions below 80% centrality. For this reason, our analysis is focusing on the data 
above 80% centrality only. 

The d+Au minimum bias collisions used in this analysis are triggered on ZDC 
detector from the east side only, which is where the Au beam is from. 

We also use the most central Au+Au collisions in the analysis which are triggered 
on ZDC detectors and the CTBs. It requires the ZDC detector with a high coincidence 
rate, and the CTBs with a large multiplicity matching the most central minimum bias 
collisions, which eventually selects the most central collisions about 12% centrality. 

2.2.4 Centrality Definition 

The collision initial condition is of great importance in heavy ion nucleus-nucleus 
collisions. The impact parameter b is defined as the distance between the geometric 
center of the two colliding nuclei in the transverse plane. The nucleons are then 
undergo interactions with each other. The number of participants referred to as N part 
is defined as the number of nucleons that participate in at least one inelastic nucleon- 
nucleon reactions. And the number of binary collisions is defined as the number 
of such inelastic nucleon-nucleon reactions, usually referred to as N co u. However, 
they cannot be measured directly from the experiment, and have to be deduced 
from other experimental observables and combined with Monte-Carlo simulations. 
The simulation is done in the geometry model with experimental measurement of 
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nucleon-nucleon cross sections and considering the multiple scattering of nucleons 
in the heavy ion collisions. Such Monte-Carlo techniques is generally referred to as 
Glauber model [27] . 
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Figure 2.5. A Glauber Monte-Carlo event (Au+Au at ^s NN = 200 GeV 
with impact parameter 6 = 6 fm) viewed in transverse plane (left) and 
along the beam axis (right). Figure is taken from |27j . 

Figure I2.5I shows a Glauber Monti-Carlo Au+Au collision event at center of mass 
v /J/v r /v = 200 GeV. The impact parameter b = 6 fm. The dark circles represent 
the participating nucleons. It is obvious that the smaller the impact parameter, the 
more participating nucleons and binary collisions in an event, hence more generated 
particles detected by the detector. Figure I2.6I shows the illustration plot of the event 
multiplicity within \r]\ < 1 distribution corresponding to the impact parameter b and 
the number of participants N part . Then experimentally, one can relate the final state 
multiplicity to the N part and b. 

In the STAR experiment, the efficiency uncorrected charged particle multiplicity in 
one unit of the pseudo-rapidity dN c h/ dr] is used instead of N c h to deduce the N part and 
b parameters, which is called reference multiplicity N Re f MuU [T21I2E]- The N Re f Mu[t 
is the number of charged particles recorded by the main TPC at mid-rapidity range 
of —0.5 < f] < 0.5. With the Nn e fMuit distribution, we can cut on certain fraction 
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Figure 2.6. An illustrated example of the total final state inclusive 
charged particle multiplicity N ch distribution with Glauber-calculated 
quantities (b, N part ). The plotted distribution and various values are il- 
lustrative and not actual measurements. Figure is taken from [27] . 



of the events to correspond to the impact parameter b and N part . Typically, we cut 
on 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70% and 70-80% from 
most central to most peripheral collisions. 
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3. DATA ANALYSIS 



In this chapter, we first give the definition of our charge multiplicity asymmetry 
variables and their correlations. We then introduce the data and quality cuts used in 
the analysis, followed by related analysis procedures and corrections. At the end, we 
check for consistency and study the systematic uncertainties. 

3.1 Charge Multiplicity Asymmetry Observables 
3.1.1 Charge Multiplicity Asymmetries 



▲ UP 




DOWN 



Figure 3.1. Schematic view of the overlap region on transverse plane in 
a typical non-central heavy ion collision. The four quadrants are labeled. 
UP = 1 + II, DOWN = III + IV, LEFT = 11 + III, RIGHT = I + IV '. 



In heavy ion collisions, the overlap area of the collision can be illustrated as an 
elliptical shape on the transverse plane as shown in figure 11.41 The event anisotropy 
can then be used to estimate the reaction-plane (RP) direction for a given event. 
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The estimation is defined as event-plane (EP). Once the event-plane is determined, 
we can then separate the overlap area into hemispheres in the collision transverse 
plane (xy plane in figure rO]) . As shown in figure l3TTj UP- and DOWN-hemispheres 
are separated by the EP. LEFT- and RIGHT-hemispheres are separated by the plane 
perpendicular to the EP. Particle multiplicity asymmetries are then defined, on event- 
by-event basis, as 

A + , UD = (N +tU - N +jD )/(N +tU + N +tD ), 
A-, UD = (N. tV - N- tD )/(N- iV + N_, D ), 
A + , LR = (N +fL - N +tR )/(N +tL + N+, R ), 

A_, LR = (AL, L - N- R )/{N- L + N- R ). (3.1) 

Here N +j u, ^+,d, and N +tR are positively charged particle multiplicities in the 
UP (quadrants I and II), DOWN (III and IV), LEFT (II and III), and RIGHT (I 
and IV) hemispheres as in figure 13.11 respectively. Those of negatively charge particle 
multiplicities are represented by N_ tU , N_^ D , N_ )L , and N_ iR . 

3.1.2 Charge Multiplicity Asymmetry Correlations 

As introduced in section [L3| the non-zero topological charges (Q w ), gauge con- 
figurations, change the quark chirality which causes the asymmetry in the number 
of left- and right-handed quarks, Nm — N R h = 2NfQ w , where Nf is the number 
of light quark flavors. Since the sign of the topological charges (Q w ) is random, 
the chirality changes are also random from domain to domain in a single event and 
from event to event. Moreover, the event-plane reconstructed from event anisotropy 
does not distinguish up from down, neither left from right. This causes the charge 
separation direction random from event to event. Therefore, only the magnitude of 
the up-down (UD) multiplicity asymmetries (A+ t uD and A- 7 ud) gets larger due to 
CME/LPV, while the sign is still random. As a result, the average asymmetries of 
UD remain zero over all events, but the distributions of A± } ud (we use A± } ud to 
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collectively denote A +t uD and A-j/d) would be wider comparing to those of A± j lr 
(i.e. ^4+,lr and ^4_,l_r), to which CME/LPV doesn't contribute. In other words, 
the variances (A 2 ± UD ) should be larger than the respective (A 2 ±LR ), where (. . .) de- 
notes the average over the event sample. Hence, we study the variance of the charge 
multiplicity asymmetries, which provide the dynamic informations about the particle 
correlations with the same charge signs. 

We also study the covariances ((A + A_)ud and (A + A_) lr) of the charge multi- 
plicity asymmetries as the opposite-sign correlations. The covariance measurement 
(A + A_) is a more traditional measurement of parity violation. The positive charged 
particle multiplicity asymmetry A + ^u£> represents, on average, the "parity- axis" of the 
orbital angular momentum and the topological charge sign; the covariance (A + A^)ud 
is then a measurement of A_ ud with respect to this "parity-axis". For local parity 
violation, the positively and negatively charged particles tend to move in opposite 
direction across the event-plane. This will generate additional negative correlation 
between the asymmetries in up-down direction comparing to left-right direction, in 
other words, (A + AJ) UD < (A + AJ) LR . 

The charge multiplicity asymmetry correlations themselves are, however, parity 
even, and subject to physics backgrounds similar to those in the charge correlator 
measurement. The backgrounds may be assessed by the left-right (LR) asymmetry 
correlations ({A 2 LR ) and (A + A_) LR ), to which CME/LPV does not contribute. This is 
because, CME/LPV is an effect along the system angular momentum direction, which 
is perpendicular to the event-plane direction. In other words, the left-right (LR) 
observables can serve as our null-reference. However, as we will show in the next 
chapter, the physics backgrounds are likely different in UD and LR measurements. 
This introduces complications in the interpretations of the measured UD and LR 
correlation differences. Details will be discussed in section 14.51 



30 



3.1.3 Dynamical Correlation and Charge Separation 

In this analysis, we compute the charge multiplicity asymmetries A±^jje> and A±^ R 
on the event-by-event basis. We obtain their variances (widths of their distributions) 
(A\ UD ) and (A\ LR ), and their covariances (A + A_) UD and (A + A_) LR . The variances 
are always positive because they are the squares of a real number. They subject to 
statistical fluctuation and detector non-uniformity acceptance effects. We have to 
subtract those effects to get the dynamical correlations, the physics we are interested 
in. 

The dynamical variances are defined as 

^(■^±,Ud) = (A±,Ud) ~ (A± t jjD,sta,t+det) > 

S(Ai LR ) = (Ai jLR ) - (Ai 

,LR,stat+det >, (3-2) 

where subscript "stat+det" stands for statistical fluctuation plus detector effects. As 
will be presented later, the dynamical asymmetry correlations of positively charged 
particles (5(A+)) and negatively charged particles (5(A_) 2 ) are consistent. We there- 
fore report the average dynamical variances (with prefix "5") as 

HA 2 UD ) = (6{Al >UD ) + 5{A 2 _ jUD ))/2, 
HA 2 LR ) = (5{Al >LR ) + 8{A\ LR )) /2. (3.3) 

We also analyze covariances (A + A-)ud and (A + A^)l R . As shown later, because 
the positively and negatively charged particles are not statistically correlated, the 
statistical fluctuation and detector effects are consistent with zero for the covariances. 
We thus do not subtract them from the covariances. 

We report the differences between the UD and LR measurements, which may be 
directly sensitive to the CME/LPV. Namely (with prefix "A"), 

A(A 2 ) ^ 5(Al >UD ) - 5(Al >LR ) « (Al >UD ) - (Al LR ), (3.4) 
A(A + A_) = (A + A_) UD - (A + A.) LR) (3.5) 
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where, ideally, the statistical fluctuation, detector effects and other backgrounds that 
are not related to the reaction plane cancel in the last step of subtraction in equation 
13.41 However, as we will show later, the detector non-uniformity causes the statistical 
fluctuation and detector effects (A\ UD stat+det ) and (A 2 ±jLRMat+det ) not identical for 
different pseudo-rapidity regions. Thus the background from UD and LR directions 
will not cancel completely. To be precise, we keep track of all the statistical fluctu- 
ations and detector effects from equation 13.21 for each variances and subtract them 
accordingly to get the dynamical correlations. And this procedure is necessary for the 
wedge size and wedge location analysis introduced in the following section, because 
the stat + det effects are largely varied with the wedge size and wedge location being 
studied, see next section. 



3.1.4 Wedge Size and Wedge Location 

We can not only investigate the charge multiplicity asymmetry between hemi- 
spheres, but also between different wedge sizes and different wedge axis locations. 



As shown in figure 13.21 we can study the charge multiplicity asymmetries with their 
correlations for any given axis with any given opening angle 2A0 W . For example in 



figure 3.2(a) , if we want to study the out-of-plane asymmetries, we count charged par- 
ticle multiplicities within the azimuthal angle relative to the EP between 90° — A0 W 
and 90° + A0 W as N±p, and that between 270° - A0 W and 270° + A0 W as N±, D . 
Similarly for in-plane asymmetries, we count particle multiplicities between 0° — A0 W 
and 0° + A0 W as N± >R , and that between 180° - A0 W and 180° + A0 W as N ±tL . Then 
the asymmetries are calculated in the same way as those defined in equation 13.11 



Furthermore, figure 3.2(b) shows the schematic configuration of the back-to-back 
wedges with axis centered at W and wedge size of 2A0 W . The asymmetries are calcu- 
lated between W ± A0 W and that between (n + <p w ) ± A0 W . In such configuration, we 
can vary the asymmetry axis W to study the progressive evolution of the asymmetries 
and their correlations from in-plane to out-of-plane. 
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We use ^4±,0 W ±A0 W to stand for the multiplicity asymmetries of positively and 
negatively charged particles with wedge axis located at W and wedge size of 2A0 W . 
Particularly, the A± >UD and A± LR we described earlier are equivalent to A±$ o± 90 o 
and ^4±,o o ±90° respectively, and we refer this configuration as hemisphere asymmetry. 




(a) Configuration for the wedge size de- (b) Configuration for the wedge location 
pendence dependence 

Figure 3.2. Schematic view of collision overlap area on transverse plane, 
and the azimuthal back-to-back wedge regions, where the charge asym- 
metries can be calculated. 

We have to be cautious when using non-hemispherical ranges. For the variances, 
it is obvious that the statistical fluctuation and detector effects do not cancel between 
the asymmetry correlations in (fi w (UD) and w+ go° (LR) wedges when we take the 
difference between UD and LR as in equation 13.41 Due to the anisotropy flow, 
we expect more particles emitted in-plane, and fewer particles emitted out-of-plane. 
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Thus, the statistical fluctuation is expected to be different for non-hemisphere wedges 
in non-central collisions, because it is approximately proportional to "1 /N" where N 
is the multiplicity within the wedge. On top of that, the electronic inefficiency 
will cause difference in the observables. The detector inefficient sector could fall into 
either UD or LR asymmetries, while in the hemisphere asymmetries which covers 
the full azimuthal angle, the inefficient sector will be guaranteed to fall into both 
UD and LR asymmetries. When calculating the UD — LR difference, the detector 
effects introduced "dynamical" correlation will not cancel in the subtraction in the 
non-hemisphere asymmetry correlations. These two factors combined together are 
non-trivial and may vary significantly for wedges with different sizes and locations. 
And this is the other reason we mentioned in last section that we want to carry 
the statistical fluctuation and detector effects, denoted as (A± i( f >w ±A<j> v ,,stat+det) > all the 
way through the analysis and subtract them respectively according to the variances. 
The statistical fluctuation and detector effects can be estimated in the same way we 
describe in section 13.71 Thus the dynamical correlation is obtained as 

w±A0 w , stat+detl ' 

(3.6) 

We report the average variance of positive and negative charged particles 

*<4w±A*,> = (*«0wW + 5« 0w± a, w » A (3.7) 

and the difference between the correlations with respect to the perpendicular asym- 
metry axises. The differences are taken between the dynamical variances after the 
subtracting statistical fluctuation and detector effects of the variances. 

For the covariances, they are not affected by such statistical fluctuation and de- 
tector effects. We thus report the raw correlations. 

3.2 Data Sets 

The data we present in this analysis were taken by STAR experiment at Brookhaven 
National Laboratory. We use minimum-bias and ZDC central triggered Au+Au colli- 
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sion data with center of mass energy y/s^N = 200 GeV per nucleon pair. We focus on 
the data taken between year 2004 and 2005 (RUN IV) with 24 million events in total. 
We also use minimum-bias triggered 200 GeV Au+Au data taken between year 2007 
and 2008 (RUN VII) for ZDC-SMD event-plane study with total 56 million events. 
As a reference, we use minimum-bias triggered d+Au data with center of mass energy 
t/snn = 200 GeV. Total events for d+Au are 9 million. The results are presented as 
a function of the number of participants, N part . 

The centrality definition is based on a Glauber Model simulation as introduced in 
section We summarize the N RefMuH cut in table O for RUN IV and RUN VII. 

Table 3.1 

Centrality definition of i/J/wv = 200 GeV Au+Au collisions and the num- 
ber of participants N part [28||29] in RUN IV and RUN VII and the average 
event-by-event < p s T<2GeV/c of RUN IV. 



# 


Centrality 




Lower N RefMu it cut (>) 


( t, 2,p T <2GeV/c) 


Npart 


RUN IV 


RUN VII 


9 


0-5% 


352.4^-4 


520 


485 


0.011 


8 


5 - 10% 


299.3+j;^ 


441 


399 


0.020 


7 


10 - 20% 


234.6+Jj 


319 


269 


0.032 


6 


20 - 30% 


166.7 + ? ° 6 


222 


178 


0.043 


5 


30 - 40% 


115.5i?i 7 2 


150 


114 


0.047 


4 


40 - 50% 


76-6^ 5 4 


96 


69 


0.044 


3 


50 - 60% 


47.825 


57 


39 


0.036 


2 


60 - 70% 


27 Atri 


31 


21 


0.025 


1 


70 - 80% 


14 1+ 3 - 6 
iq: - i -5.0 


14 


10 


0.017 
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3.3 Quality Cuts 

In order to get quality event and track information, we apply STAR standard 
cuts to ensure the events and tracks are with good precision, yet least biased by the 
detector imperfection. The cuts used here are described in previous chapter. 

3.3.1 Event Selection 

Event wise, a minimum-bias trigger (section 12 . 2 .3[) was used at data taking. Fur- 
thermore, events are required to have the collision vertices within 30 cm from the 
STAR detector center along the beam line (\v z \ < 30 cm), in order to ensure the 
events are not biased toward one side of the TPC. 

To reject collisions from beam halo and beam pipe, or collisions from possible 
secondary vertices, we cut on the radius of the vertex position to the beam pipe 
center with v r = ■ s /v% + Vy < 2 cm. 

During data taking, the TPC magnet operated in two polarity configurations, full 
field (FF) and reverse full field (RFF), with magnetic field strength of 0.5 T in beam 
line direction. We then use different magnetic polarity to assess the detector effects 
which is effectively the same as switching positive and negative charges. 

The reference multiplicity of an event is also required to be less than 1000 in order 
to reject pile-up events. 

3.3.2 Track Selection 

Track wise, tracks are required to be within ±1 in pseudo rapidity (\rj\ < 1), 
where the TPC has the best tracking performance. Each track is required at least 20 
hit points (out of 45 at most) used in track reconstruction from the TPC. Also, the 
ratio of the number of hit points in track reconstruction to the most possible number 
of hit points is required to be larger than 51%, which will eliminate partial track 
reconstructed from a single track. We also require a lower transverse momentum cut 
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with pt > 0.15 GeV/c, which is the lower limit of STAR detector. We require each 
track with DCA < 2 cm to ensure the particle is from the collision vertex. 

We vary the event cuts and track cuts to study the systematic uncertainties, which 
is discussed in section 13.91 

We summarize the cuts in the table 13.21 and the number of events after the cuts 
for each dataset in table 13.31 



Table 3.2 

The standard event and track selection cut. 



Event selection cut 



Vertex z position v z \v z \ < 30 cm 

Vertex radius cut v r = x/v% + v% v r < 2 cm 

Reference multiplicity N RefMu it N RefMu it < 1000 
Track selection cut 

Pseudo-rapidity t] \rj\ < 1 

Number of hits nfit nfit > 20 

Ratio of nfit to maximum fit points rfit > 0.51 

Lower transverse momentum p^ > 0.15 GeV/c 

DCA DCA < 2 cm 



Table 3.3 

Dataset and statistics. 

Dataset Number of events after cuts 

RUN IV Au+Au 200 GeV min-bias 22.5M 

RUN IV Au+Au 200 GeV top 2% central 5.5M 

RUN IV d+Au 200 GeV min-bias 8.9M 

RUN VII Au+Au 200 GeV min-bias 56.4M 
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3.4 Detector Efficiency Correction 

STAR detector has full azimuthal coverage in 2tt. For a perfect detector, the 
particle distribution accumulating over a large amount of events should be uniform 
with respect to the azimuthal angle. But in reality, the TPC is made of 12 separated 
sectors on each side of the endcaps. The sector boundaries and deficit sectors due to 
electronics failure will reduce tracking efficiency, and create "dynamical" correlation 
on the event-by-event basis. To reduce the "dynamical" correlation introduced by 
non perfect detector, we will have to correct the single track efficiency to flatten the 
azimuthal angle distribution. 



Figure 3.3(a) shows an example of the lab frame azimuthal angle distributions 
of positively charged particles before the acceptance correction. Data used in the 
figure are from RUN IV Au+Au 200 GeV collisions from east-side (77 < 0) in red 
and west-side TPC (77 > 0) in blue separately. The centrality range is 30-40%; for 
all centralities and charges, refer to figure IA.1I The magnetic field polarities and the 
different charges have been summed together and the transverse momentum px is 
integrated over 0.15 < pt < 2 GeV/c. From the uncorrelated raw distributions, the 
repeatedly dropping pattern, especially in most central collisions with large statistics 
flA.lj) . is clearly seen because of the sector boundaries of the TPC. And there is 
additional inefficiency in east-side TPC (r] < 0, blue), range within 5ir/3 < <fi < 2-k. 
This is due to the inefficiency in the electronics readout system (RDO) of two out of 
twelve sectors in the east-side TPC. We found the effect persists over the entire data 
set of RUN IV period, with no significant time variation. 

We correct for the single particle dependent inefficiency, mainly due to sector 
boundaries and electronic dead sectors, by giving each particle a weight depending on 



the particle position. We normalize the raw <fi distribution in figure 3.3(a) to average 
unity. The normalized <ft distributions are defined as acceptance x efficiency, e(4>). 
Those e(<fr) distributions are separated according to magnetic field polarities, particle 
charges and centrality bins. They are then further separated for positive and nega- 
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Figure 3.3. Example plot of detector acceptance effect for positively 
charged particles in 30-40% Au+Au 200 GeV RUN IV collisions. Panel 
(a) shows the raw (uncorrected) lab frame azimuthal angle (0) distribu- 
tion. Blue curve is for those particles in rj < (east-side TPC) region, and 
red curve is for particles in > (west-side TPC) region. Data include 
both full field and reverse full field magnetic polarities, and the transverse 
momentum is integrated over 0.15 < px < 2 GeV/c. Panel (b) shows the 
corrected distribution from the same data sample. 
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tive r) (corresponding to east- and west-side of the TPC tracks), positive and negative 
vertex in z direction (v z ), and for different bins as following, 0.15 — 0.5, 0.5 — 1.0, 
1.0 — 1.5, 1.5 — 2.0 GeV/c. Then, we take l/e(0) as the single particle weight as our 
correction factor for the detector acceptance x efficiency. Note, the weight l/e((f>) 



could be either larger or smaller than unity. Figure 3.3(b) shows the corrected single 



particle </> distributions of the same set of particles in figure 3.3(a) The corrected 
distributions are almost uniform, although some jitter effects are still not completely 
removed due to binning and fluctuation issues, which means we compensate the in- 
efficiency in the single particle <fi non-uniformity. The weight l/e(0) is then applied 
to all the asymmetry multiplicities and the TPC event-plane reconstruction through 
this analysis. For all centralities and charges, the corrected <fi distributions are shown 
in appendix figure IA.2I 

There are two main reasons we have to apply the single particle <p weighted correc- 
tion. The first reason is to flatten the event-plane reconstructed from the TPC tracks. 
By nature, event-plane direction is random and should be flat in azimuthal direction. 
As we will introduce in section 13.5. 1\ the second order event-plane reconstruction 
method is based on particle azimuthal distribution anisotropy. If we do not correct 
the detector non-uniform </> inefficiency, the artifact anisotropy will couple with the 
particle distribution anisotropy making the reconstructed event-plane have preferred 
azimuth direction. We do see this effect in data for about several percent magnitude. 
After we applied the <fi weighted correction, the event-plane distribution is flat for all 
centralities. We will show that in section 13.5.11 

The other reason for dependent efficiency correction is to compensate the bias 
of charge asymmetry correlations due to the detector effect. Figure 13.41 shows the 



selected charge asymmetry correlations in UD examples: (A +UD ) (3.4(a)), (A_ UD ) 



(3.4(b)) and (A+A_) (3.4(c)), before and after the dependent efficiency correction. 
In the plots, all the asymmetry correlations are multiplied by the number of partic- 
ipants Np ar t to better show the magnitude. Note, asymmetries are calculated with 
tracks from half side of the TPC 77 region (77 > or 77 < 0) with respect to the EP 
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Figure 3.4. Asymmetry correlations: panel (a) {A\ UD ), (b) (A 2 _ UD ), 
(c) (A + A-)ud (scaled by number of participants, N part ), before and after 
single particle corrections for the <fi dependent acceptance x efficiency 
The EP is reconstructed by charged particles with p T range of 0.15 < px < 
2.0 GeV/c from one side of the TPC, and the asymmetry correlations are 
calculate in the same pt range but from the other side of the TPC. The 
lines are the statistical fluctuation and detector effect estimation. For 
clarity, only selected asymmetry correlations are shown. For all plots, see 
appendix. 



reconstructed with tracks from the other half side of the TPC r\ region (77 < and 
77 > 0). The reason for this setup is to avoid self-correlation, which is detailed in 
section 13.61 

We show only selected charge and UD correlations in the figures for clarity reason. 
The efficiency correction effect for other charge asymmetry combinations and corre- 
lations in LR is similar to UD. They are shown in figure IAT41 and IA~5l We see greater 
effect on correlations for asymmetries from rj < region than that from r\ > region. 
This is because of the electronics RDO system inefficiency within 5tt/3 < <fi < 2tt 
for east-side of the TPC (77 < region). The inefficiency creates additional non- 



41 



uniformity, and more "dynamical" correlations from the east-side than the west-side 



of the TPC by comparing the difference before correction in figure 3.4(b) to 3.4(a^ 



The correction effect is larger in the opposite-sign correlations ((A + A_) in figure 
3.4(c)) rather than the same sign correlations {(A 2 ) in figure 3.4(b)). As seen in 



figure 3.4(c), the {A + A_)jjd for the r\ < region is significantly larger than the 
i] > region before the <fi efficiency correction. After the single particle efficiency 
correction, the correlations are consistent between different charges, and different r\ 
regions, which are shown in section 13.81 The lines in the plots will be discussed in 
section 13.71 



3.5 Event Plane Reconstruction 

In this section, we introduce two methods to reconstruct the event-plane. The 
second order event-plane is reconstructed from the TPC tracks. And the first order 
event-plane is reconstructed from the ZDC-SMD neutron energy shower profiles. We 
will also discuss the event-plane resolution, which defines how well the estimation is 
for the reconstructed event-plane ipEP compare to the real reaction plane ippp, even 
we don't know it exactly. 



3.5.1 Second Order Event Plane Reconstruction from TPC 



As we introduced in section 11.41 the initial geometry eccentricity of a collision 
event will translate into the final state particle azimuth and transverse momentum 
distribution anisotropy, which is called flow. The anisotropic flow can be described 
by a series of non-zero Fourier coefficients shown in equation 11.51 Then the reaction 
plane can be estimated from each term of the Fourier components. 

In general, for a given event, its n-th order reaction-plane can be estimated in the 
following way [3"U] : 

Y Wi sin(n0, 

ip n ,EP = I arctan — — I /n. (3.8) 
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The sums go through all the particles used in the reaction-plane reconstruction. Wi 
is the weight used to get the optimized reaction plane resolution for each particle. 

In this analysis, we mainly focus on the study of asymmetry correlations with 
respect to the second order event-plane, which gives the best reaction-plane ip2,RP 
estimation, because the second order term in the Fourier expansion is the dominant 
term of the event shape anisotropy. For better event-plane resolution, the particle 
transverse momentum px is used as the weight, Wi = pt,i because t>2 increases with 
Pt [31]. Also we have to apply the dependent single particle efficiency correction. 
The second order event-plane i/)ep is then calculated as 

Y,Pr,i sin(20i)/e 
* B „-|«c*»£ |/2. 
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Figure 3.5. Reconstructed second order event-plane azimuthal distribu- 
tion for RUN IV Au+Au 200 GeV collisions in 30-40% centrality. The 
event-plane is reconstructed from charged particles within 0.15 < pr < 2 
GeV/c from 77 < (red) and rj > (blue) separately. The lines are event- 
planes reconstructed with particles without single particle <fi acceptance 
correction, and data points are with correction. 
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Figure [3~5l shows the reconstructed event-plane azimuthal angle ipEP distributions 
example from RUN IV 200 GeV Au+Au collision in 30-40%. The complete event- 
plane distributions of all centralities are shown in figure IA.3I 

Before applying the single particle <j) acceptance correction, the event-plane dis- 
tributions are shown as lines. The distribution of event-plane reconstructed with the 
particles from east-side TPC within electronic deficit region (rj < 0), shown in dashed 
blue line, has a larger deviation from flatness compared to the other side. This is 
because the EP reconstruction method based on the particle distribution anisotropy 
will try to find a direction with more particles emitted in the event-plane direction, 



and fewer particles emitted in out-of-plane direction. And as shown in figure 3.3(a) 
the inefficient sector in r\ < region is located at 57r/3 < <ft < 2n in azimuth angle. 
The inefficiency gives a preferred direction to the event-plane at approximately 7r/3, 
which is perpendicular to the center of the deficit sector around ll7r/6. The peak 
of the blue line in figure [3751 is found at 7r/4, which is roughly consistent with our 
expectation. 

To correct such detector effects, we apply the single particle <fi dependent correction 
to all the tracks used in the EP reconstruction. The results are shown as data points 



in figure I375j which are flat for both 77 > and r\ < regions. So, we are confident with 
the reconstruction method and the necessary correction applied to remove detector 
non-perfection effect in the EP reconstruction. 

Although we corrected the detector effects in event-plane reconstruction, the re- 



constructed EP is still an approximation of the true reaction-plane. In equation I3.8[ 
the event-plane will approach to the true reaction-plane if the event had infinite num- 
ber of particles. However, a real event has only finite number of particles recorded in 
TPC. This will limit the accuracy of the estimated reaction-plane. The inaccuracy is 
defined as event-plane resolution [30] 



e EP = (cos(2(^ £ p - ip R p))), (3.10) 
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for second order event-plane ipEP and reaction plane i/jrp here. Under such definition, 
€ep is 1 if the ipEP is exactly the same as ippp, and if ipEP is completely random 
relative to ippp. 

Note we separate an event into two sub sets, particles from east- and west-side of 
the TPC tracks. By nature, the RP is a collision parameter which should be identical 
for both sub sets because they are from the same event. However, the event-plane 
reconstructed from the two sub sets may not be the same. And it is even possible 
that, due to other physics mechanism such as resonance decay or jet quenching effects, 
there are events with odd shapes making the event-planes significantly different from 
one another. Although we do not know the true value of ipRPi t ne event-planes 
reconstructed from the sub events are subject to the same resolution effects with 
respect to the true reaction-plane. We can thus show that the event-plane resolution 
can be calculated from the sub events approximately as [30J , 

(COS 2 (tpEP, v >0 - 1pEP,r,<o)) = (COS (2 {lpEP, v >0 ~ ^Rp) ~ 2 (lpEP, v <0 ~ IpRp))) 



where ipEP,r]>o and ipEP,rt<o are the reconstructed event-plane azimuthal angles from 
particles in 77 > and 77 < regions respectively. Note, this particular event-plane 
resolution comes from the sub-events from the same event, which is most relevant for 
the studies we carry out in this thesis. 

Figure 13.61 shows the EP resolution €ep as a function of centrality (number of 
participants N part ). The maximum event-plane resolution is found in medium central 
collisions. The resolution decreases in the peripheral and most central collisions. This 
is because the event-plane reconstruction method is essentially based on the particle 
azimuthal anisotropy which defines the event-plane separating the particles into equal 
halves. In very peripheral collisions, the multiplicity is low. Those events are more 
likely being affected by non-flow and fluctuation. For example, a pair of back-to-back 



« (cos 2 {ip E p,r,>a - i>Rp) cos 2 {ipEP,r,<a - *Prp)) 



e EPi 




(3.11) 
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Figure 3.6. Second order event-plane resolution from sub-events (77 < 
and rj > 0) as a function of the number of participants in RUN IV Au+Au 
200 GeV collisions. Statistical error bars are smaller than the symbols. 



di-jets may define the event-plane direction, which is unrelated to the reaction-plane. 
Thus the resolution is low. On the other hand in the most central collisions, the two 
colliding nuclei have the maximum overlapping area, less anisotropy in other words. 
Thus the EP resolution is also low. In the medium central collisions, they have more 
significant anisotropy, resulting in the largest EP resolution. 

After reconstructing the event-plane, correcting for the detector effects and having 
the EP resolution under control, we want to study how the EP resolution affects our 
charge asymmetry variances and covariances. We cannot reach perfect EP resolution, 
but we can reduce the resolution by randomly throwing away a certain fraction of 
the particles during the event-plane reconstruction. In this way, we can show and 
understand how the asymmetry correlations vary with different EP resolutions. 

Figure 13.71 shows EP dependences of the 20-40% mid-central Au+Au 200 GeV 
RUN IV collisions. The rightmost data points in each figure are asymmetry correla- 
tions with respect to the event-plane reconstructed using all charged particles from 
the half side of the TPC. Then from right to the left, we artificially reduce the EP res- 
olution by discarding 25%, 50%, and 75% tracks used in the EP reconstruction. The 
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EP resolution is estimated by e^p = a/ (cos 2 (ipEP,r]>o — ipEP,r)<o))- We then plot the 



charge asymmetry variances (A 2 ) (figure 3.7(a)), covariances (A + AJ) (figure 3.7(b)) 



and their differences between UD and LR directions (figure 3.7(c)) as a function 



of the corresponding EP resolution. Note that the asymmetries are calculated with 
all charged particles from half side of the TPC. We do not discard particles in the 
asymmetry calculation. 
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Figure 3.7. Charge multiplicity asymmetry correlations: (a) (A 2 ), (b) 
(A + AJ) and (c) their differences between UD and LR, as a function of 
the event-plane resolution €ep in 20-40% mid-central centrality. The solid 
lines are linear fits to the data. The dashed lines are linear fits with fixed 
zero intercept at €ep = 0. Error bars are statistical. 
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In figure 13. 7\ we apply a linear fit of the correlation data to extrapolate the EP 
resolution to zero and unity. The fits are presented as solid lines for all the variances, 
covariances and the UD — LR differences. When EP resolution approaches zero, i.e. 
the event-plane is random, we cannot distinguish between UD and LR. Then the 
asymmetry correlations of UD and LR should converge at the same intercept. This is 



shown in figure 3.7(a) and 3.7(b) , where the linear fits of variances and covariances in 
U D and LR roughly converge at the zero EP resolution. The differences between the 
UD and LR, A(A 2 ) and A(A + A_) correlations should vanish at zero EP resolution 



because UD and LR are identical at zero EP resolution. This is shown in figure 3.7(c^ 



where the linear fit in solid lines roughly converge in zero at zero EP resolution. The 
dashed lines in the same figure are the fits with fixed intercept of at zero EP 
resolution, which are visually consistent with what we expected. 

The linear fits in all three figures indicate that EP resolution only smears out the 
correlations between UD and LR. If we extrapolate the EP resolution to unity, the 
magnitude of the correlation differences between U D and LR will only be larger than 
what we actually measured from data. In this analysis, we do not correct for EP 
resolution. This is because, the linear extrapolation only works when the high-order 
harmonic terms of asymmetry correlations in equation 14.51 are negligible, which is not 
true for the variances as we shall discuss in section 14.5.11 The high-order terms of the 
variance contribute significantly to the UD — LR correlation. It remains unknown 
how the high-order harmonic terms respond to the EP resolution. Similar figure IA.6I 
shows the most central and peripheral charge asymmetry correlations and their UD 
and LR differences as a function of the EP resolution. And more detailed plots of each 
centrality are shown in figure IA.7I for variances, figure IA.8I for covariances and figure 
IA.9I for the differences of UD and LR. It is important to realize that our qualitative 
conclusions will not change if we have perfect EP resolution. 

Note that, the second order event-plane azimuthal angle ipEP range from to 
Ti. The event-plane angle ipEP is equivalent to ipEP + 7r. When calculating the 
asymmetries, we randomly flip the reconstructed event-plane to make it range from 
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to 2tt. This is because, if we do not flip ipEP, the event-plane will have a preferred 
azimuthal direction in one half side of the TPC. Thus the UD and LR will also 
have a preferred direction in azimuth. This will introduce systematic errors due to 
any residual effect from the imperfect detector efficiency. After random flipping, the 
preferred direction is avoided, so is the systematic uncertainty due to the preferred 
direction. 

3.5.2 First Order Event Plane from ZDC-SMD 

The Zero Degree Calorimeters and Shower Maximum Detectors (ZDC-SMD) are 
located ±18 m away from the center of the STAR detector. It records the neutron 
energy deposit profile from the deflected spectators, which can be used to measure the 
first order event-plane determined by direct flow. Since its pseudo-rapidity coverage 
(Ir/I > 6) is far away from TPC (|t/| < 1), and the measured neutrons are from 
the fragmented gold nuclei which do not participate in the reaction, there is little 
correlation between the ZDC-SMD signal to TPC tracks. By using the event-plane 
reconstructed from ZDC-SMD, we can further remove possible physics correlations 
between the charge asymmetry correlations and the event-plane. However, the first 
order event-plane resolution is not as good as the TPC event-plane as we will show 
below. 

The ZDC-SMD detector measures the energy deposit profile with 7-slate vertical 
and 8-slate horizontal channels from both east-side and west-side of the STAR de- 
tectors. The raw signals are corrected by pedestal subtraction and electronic gain 
corrections for each channel. The pedestal subtraction is applied to correct the differ- 
ent electronic background of each readout channel. And the gain correction is used 
to correct the linearity of the ADC response to the neutron energy deposit. After 
the corrections, the signals give a hit profile in the transverse plane in the manner 
of vertical and horizontal distributions. The vector from the beam center (averaged 
for each run) to the profile center gives the direct flow direction of the collision on 
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each side of ZDC-SMD in transverse plane, which is the measurement of the first 
order event-plane direction from one side of the ZDC-SMD. Combining the two vec- 
tors from east- and west-side, we can get a better measurement of a single first order 
event-plane for each event [32J. Note that the two vectors are mostly back-to-back 
due to momentum conservation of the fragmented spectators. Also note that, the 
first order event-plane ranges from to 2tt, while second order event-plane ranges 
from to 7r. 
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Figure 3.8. Reconstructed ZDC-SMD first order event-plane distribu- 
tion for RUN VII Au+Au 200 GeV 30-40% centrality. The blue curve 
shows the uncorrected raw EP distribution. The red curve shows the cor- 
rected EP distribution with the recentering method. The event-plane is 
reconstructed from combining the shower profile from east- and west-side 
ZDC-SMD. Error bars are statistical only. 



Figure [3781 blue curve shows the raw ZDC-SMD first order event-plane of RUN VII 
200 GeV Au+Au collision at 30-40% centrality. As we can see, the raw event-plane 
distribution is largely non-uniform in the azimuth, and has the preference direction 
with a fluctuation as large as about 50%. Thus the ZDC-SMD event-plane has to be 
corrected for non-uniformity in the azimuthal direction same as for the second order 
event-plane. The major difference between these two is that, the TPC event-plane 
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reconstruction is track based, while ZDC-SMD event-plane reconstruction is profile 
based. We can correct each track to flatten the TPC event-plane, but it is impossible 
to do the same to the energy profiles. So another method, the so-called "recentering" 
method, has to be used for the analysis [33] . 

The idea is to shift the event-plane angle i^ep by a correction according to its 
location 5iJj(iJjep), such that the new event-plane angle distribution ip' EP = ip E p + 
8ip{ipEp) is flat for the whole event sample. The raw event-plane distribution can be 
Fourier decomposed as the following: 

dN(i/j E p) a 



dip 



EP 



+ ^2 [ a ™ cos (?#£:p) + b n sm(nip E p)}, n > 1, 



with 



7i 



a n = - N(ip E p) cos(mp E p) dip E p, n > 0, 



o 



1 f 2n 

b n = - N(ip E p) sin(nx) dip E p, n>\. 
Jo 



In order to obtain a flat event-plane distribution, the correction term must satisfy the 
following: 

1p' EP = IpEP + Slp(lp E p), 

dAT(W) = oo 
d^P ' 2 ' 

which leads to 

8ip(ipEp) = 2^ cos ( n ^£p) + B n sin(nij; EP )}, n > 1. 

n 

We can then easily derive that 

2 

A n = (sm(nij; E p)), 

n 

2 

B n = -(cos(mp E p))- 
n 

In principle, if we apply infinite number of orders (n) to the correction, we defi- 
nitely will get a flat first order event-plane. In this analysis, we take the shifting up 
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to the 4th order correction (n = 1, 2, 3, 4), and the corrected event-plane distribution 
is shown in figure 13.81 in red data points. With up to the 4th order correction, the 
final EP distribution is flat enough in the azimuthal angle. Also in figure lA.lOt 
we show all centrality first order event-plane azimuth distributions before and after 
recentering correction. The data we used are RUN VII 200 GeV Au+Au collisions. 
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Figure 3.9. The first order event-plane resolution as a function of the 
number of participants in RUN VII Au+Au 200 GeV collisions. Statistical 
error bars are too small to be seen. 

To show how good the first order event-plane is compare to the reaction-plane, we 
calculate the event-plane resolution which is similar to the second order event-plane. 
It is calculated as: 

C£P = \j^{cOs{lp E p^DCEast ~ ^EP,ZDCWest)) ■ {3.12) 

Note that we combine the two event-planes 4>EP,ZDCEast and ipEP,ZDCWest to form a 
single event-plane, so the resolution has a y/2 difference to the second event-plane. 
Figure 13.91 shows the first order event-plane resolution as a function of N part for 
Au+Au 200 GeV RUN VII data. It is lower than the second order event-plane shown 
in figure 13.61 
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3.6 Self-Correlation 

In general, when working on correlation study, one has to be cautious about self- 
correlation. The problem rises when we calculate an observable from one set of data, 
and then correlate it with another observable calculated from the same set of data. 
The two observables are intrinsically related, and the correlation between these two 
is automatically affected by self-correlation. 

Particularly in this analysis, the reconstructed second order event-plane method 
utilizes the particle distribution anisotropy. The EP reconstruction algorithm guar- 
antees that the EP always divides the event multiplicity into more or less two equally 
halves in UP and DOWN hemispheres. For example, in one set of particles from 
the same event, if the positively charged particle multiplicity is unbalanced toward 
one side of the event-plane either due to fluctuation or underlying physics, the neg- 
atively charged particle multiplicity is more likely unbalanced toward the other side 
of the event-plane. Therefore, the asymmetries of positively and negatively charged 
particles are anti-correlated between UD with respect to the EP reconstructed from 
the same set of particles. However it does not affect the correlations between LR. 
This results in smaller covariance {A + AJ)jj£, than (A + A_) m, and the difference is 
an artifact of self-correlation, which does not indicate physics dynamics. 

To show the self-correlation effect, we calculate four sets of asymmetry correlations 
and compare them. They are separated according to rj regions for EP reconstruction 
and asymmetry calculation in different combinations. 

(I) Using particles within —1 < rj < for EP reconstruction and within < r\ < 1 
for asymmetry correlations; 

(II) Using particles within < 77 < 1 for EP reconstruction and within — 1 < 77 < 
for asymmetry correlations; 

(III) Using particles within — 1 < 77 < for EP reconstruction and within — 1 < r] < 
for asymmetry correlations; 
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(IV) Using particles within < rj < 1 for EP reconstruction and within < rj < 1 
for asymmetry correlations; 

The results of the four cases are shown in figure 13.101 for the covariances and figure 
13.111 for the variances. 
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Figure 3.10. Panel (a): The covariances of the four cases in UD direc- 
tion. Panel (b): The covariances of the four cases in LR direction. The 
covariances are scaled by the number of participants N part and plot against 
centrality (N part ) for four cases. The particle px range of 0.15 < pr < 2.0 
GeV/c is used for both EP reconstruction and asymmetry calculation. 
Error bars are statistical. 



To better illustrate the differences, figure 13.121 shows the relative ratios between 



UD and LR for the covariances {A + AJ)ud/ (A + AJ) lr in figure 3.12(a) and the vari- 



ances (Ajj D ) / (A 2 LR ) in figure 3.12(b) with the above four combinations. The signifi- 



cant self-correlation effect is observed in the covariances from figure 3.12(a) of cases 
(III) and (IV), where the ratios of the cases with the asymmetry correlations and 
event-plane reconstructed from the same side of the TPC are diverged from those 
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Figure 3.11. Panel (a): The variances of the four cases in UD direction. 
Panel (a): The variances of the four cases in LR direction. The variances 
are scaled by the number of participants N part and plot against centrality 
(N part ) for four cases. The particle pt range of 0.15 < px < 2.0 GeV/c is 
used for both EP reconstruction and asymmetry calculation. Error bars 
are statistical. 



cases using different side of TPC tracks. The effect is more significant for the pe- 
ripheral events, where the total multiplicity is less than in the central events. Thus, 
the fluctuation to multiplicity ratio due to self-correlation is larger in the peripheral 
events than in the more central events. However, the ratios from different side of 
the TPC tracks are relatively stable over all centralities, which indicates a smaller 
self-correlation effect. In this analysis, we use the average of cases (I) and (II) as our 
result correlations to avoid the self-correlation. 

On the other hand, the same-sign correlations, the variances (A 2 ± UD ) and (A\ LR ) 



show no significant self-correlation effect as shown in 3.12(b) This is because the 
particles used in variances and the EP reconstruction are not identical set of particles. 
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Figure 3.12. Relative ratios of charge multiplicity correlations as a 
function of the number of participants N part for four combinations of rj 
ranges used in EP reconstruction and asymmetry calculation. Panel (a): 
(A + A_) UD /{A + A_) LR , Panel (b): (A 2 UD )/(A 2 LR ). The particle p T range of 
0.15 < pt < 2.0 GeV/c is used for both EP reconstruction and asymmetry 
calculation. Error bars are statistical. 



Only half of the particles (positive or negative charged particles) are used in the 
asymmetry calculation. Thus the self-correlation effect is not obvious from the ratio 
plot. We still use the average of cases (I) and (II) for the variances to consistent with 
the covariances. 

In such setup, we calculate the charge asymmetries and their correlations from 
one side of the TPC tracks (from 77 > or 7] < 0) with respect to the event-plane 
reconstructed from the other side of the TPC tracks (from 77 < or 7] > 0). For 
each event, we have two sets of asymmetries and their correlations, and we take the 
average of them to increase statistics. 
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To further remove the short range correlation which has a bulk correlation with 
the rapidity span around one unit for the soft particles, the ZDC-SMD first order 
event-plane is used. It extends the pseudo-rapidity range to \t]\ > 6. The correlation 
between the TPC tracks and the ZDC-SMD signal are reduced to minimum. In order 
to make direct comparison, we still divide an event into two sub-events according to 
the pseudo-rapidity range for the asymmetry and their correlation calculation with 
respect to the full event-plane reconstructed from ZDC-SMD. 



3.7 Statistical Fluctuation and Detector Effect 

It is obvious that the variances are non-zero even if there are no dynamical fluc- 
tuation presents, because they are defined as the squares of the ratios which are real 
numbers. This is a trivial effect of the statistical fluctuation of the finite multiplic- 
ity. Take N± t u = (N± t u) + 5N± } u as the collective notation for N +t u and N- t u, and 
N±.d = (N± t D} + 5N±,d for N +t £> and AL^. If we assume their fluctuations are Pois- 
son, the statistical fluctuation effect of the variances can be expanded in 5N±/N± 
and approximated by 

2^ 



(A±UD,stat) 



5N±n - 6N±n 



(NQ \ + {6N ±>U + 6N ±!D ) /{N±) , 
(N±) + 1 



(3.13) 



where (N±) = (N±p) + (N±,d), and likewise for (A 2 ±LRstat ). The multiplicities 
in equation 13.131 are the measured multiplicities before any acceptance corrections. 
The acceptance corrections cancel in the correlations and do not contribute to the 
statistical fluctuations. Example results are shown for (A 2 + LR ) in rj < region in 



figure 3.13(a) and rj > region in 3.13(b) as black curves. The method is referred as 



"1/iV" approximation. The data are scaled by the number of participants N part to 
better show the difference. The approximation works when the multiplicity is large 
enough, and the multiplicity distribution is close to Poisson distribution. 
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Figure 3.13. Panel (a): Statistical fluctuation and detector effects in 
charge asymmetry variances scaled by the number of participants N part 
from east-side of the TPC, rj < region, with respect to the EP recon- 
structed from west-side of the TPC, 77 > region. The black curve shows 
the "1/iV" approximation by equation 13.131 The blue curve shows the 
pure statistical fluctuation (A 2 + LR stat ) with "50-50" method. The statis- 
tical fluctuation plus detector effects (A? + LR . stat+det ) are shown in green 
crosses using scramble method, and red circles using flipping- tt method. 
Panel (b): Same as (a) but for rj > region with EP reconstructed from 
rj < region. Panel (c): The ratios of statistical fluctuation estimations 
between different methods of different rj regions. Panel (d): The consis- 
tency check of ratios of UD to LR and positive to negative charge obtained 
by adding- 7r method. Data are from RUN IV Au+Au 200 GeV collisions. 
The particle pt range is integrated over 0.15 < pr < 2.0 GeV/c. 
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The Poisson assumption of the multiplicities may not be true in real data. To 
improve the estimate, we can count the multiplicities of each charges, and then obtain 
the statistical fluctuation from real data. It is done by generating a random number 
with 50% chance to flip a particle to its opposite direction in azimuthal angle, in 
other words, randomly flipping a particle to up or down hemisphere, and left or right 
hemisphere, without considering single particle efficiency correction l/e(0). After 
the random flipping, the asymmetries and their correlations are calculated again as 
the statistical fluctuation and noted as (A 2 ±stat ). Under such process, the charge 
correlations are completely destroyed within (A± stat ) because the flipping is random 



and does not depend on charges. The result has been shown in figure 3.4(a) and 



3.4(b) without 0-correction. We refer this method as the "50-50" method, which 



doesn't require the multiplicities to be Poisson distributions or very large numbers. 
After the random flipping, the single particle distribution does not reflect the detector 
acceptance efficiency; in other words, it would work only when the detector is perfectly 
uniform, because the flipping does not dependent on the detector performance. Thus, 
this method is an estimation of the pure statistical fluctuation only. The blue curves 



in figure 3.13(a) (rj < 0) and 3.13(b) (rj > 0) are the statistical fluctuations obtained 
by this method from the same sets of data used in the "1/iV" method. As shown in 
the figures, the "1/iV" method underestimates the statistical fluctuation especially 
in low multiplicity events, because the large multiplicity and Poisson distribution 
assumptions cannot describe the data precisely in those events. 

The ratio of the statistical fluctuations between "50-50" method and U l/N" method 
is shown as the red line from rj < region and blue line from rj > region on fig- 
ure 3.13(c) to quantify the difference between the two methods. The "1/JV" method 



underestimates the statistical fluctuation for about several percent in peripheral to 
mid central collisions, and agrees with "50-50" method very well in the most cen- 
tral collisions. It is because the multiplicity Poisson distribution assumption breaks 
down in low multiplicity events, where the "1/iV" method underestimates. In high 
multiplicity events, the U l/N" method works well where the Poisson distribution 
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holds. The rj > and r/ < regions show similar effect of the multiplicity off Poisson 
distributions. 

The "50-50" approximation assumes a perfect detector with uniform distribution 
in the azimuthal direction. However, our detector performance is subject to the 
sector boundaries, hardware efficiencies and other deficit electronics, etc. As men- 
tioned in section 13.41 the non- uniformity has to be corrected for the asymmetries and 
event-plane reconstructions with the average <fi dependent correction factor l/e(0) 
respectively. The correction factor can correct the single particle distribution and 
the event-plane distribution, however, the event-by-event "dynamical" fluctuations 
introduced by the asymmetry of detector inefficiencies cannot be corrected. We refer 
detector induced event-by-event "dynamical" fluctuations as "detector effects" , which 
have to be corrected to obtain the physics dynamics. To assess the detector effects, 
we modify the "50-50" method to account for the detector non-uniformity. 

To do so is, for each particle, we still flip its azimuthal angle by adding n to its 
azimuthal angle, or do nothing. But instead of 50% chance of flipping as we did in 
"50-50" method, the probability for adding ir to the azimuthal angle is determined 
by the relative azimuthal acceptance x efficiency l/e(4>) and l/e(0 + 7r). For example 
a particle with azimuthal angle <ft, the chance of flipping it to its opposite direction is 
e(0 + 7r)/(e(0) + e((j) + ir)). We then calculate the asymmetries and their correlations 
after the random flipping and note that as (A 2 ± stat+c i et )- The result has been shown in 



figure 3.4(a) and figure 3.4(b) with 0-correction. The method is referred as adding- tt 



method, and the selected correlations are shown with red circles in figure 3.13(a) for 



{A 2 + lr stat+det) of ?7 < region and figure 3.13(b) of r\ > region. The magnitude of 
(A^tat+det) i s larger than the pure statistical fluctuation (A 2 stat ). We show the ratios 



of the two in figure 3.13(c) with black (77 > region) and green (77 < 0) lines. In 
the figure, we can see that the detector effect is a few percentage larger than the 
pure statistical fluctuation. The detector effect is larger in the most central collisions 
because the particle multiplicity is larger in more central collisions, therefore, the 
detector effect and tracking inefficiency are more significant. Also note, the ratio of 
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7] < region in green line is slightly larger than the black link which is the r/ > region, 
because the detector deficit electronics effect induces the "dynamic" correlations in 
this region. 

Using Monte-Carlo simulation is another way to assess the statistical fluctuation 
plus the detector effects. For a given event, we generate the same amount of particles 
as the TPC recorded with random azimuthal angles according the corresponding 
distribution. During the particle generation, we keep other parameters unchanged, 
such as pt and charge. We use the generated events as mixed event to calculate 
the charge asymmetries and their correlations. Since the Monte-Carlo simulation 
does not depend on the particle charges, the mixed event asymmetry correlations 
are actually the measurement of the statistical fluctuation. Because we also use the 
detector acceptance to mimic the detector effect, the asymmetry correlations from the 
mixed events also take care of the detector effects. We refer this method as scramble 



method, and the selected results are shown in fig ure |3. 13(a)] for {A^_ stat+det) ^ 



rj < region and in figure 3.13(b) in 77 > region with green crosses. The results 
show good consistency with the adding-7r method for all charge, 77, and UD and LR 
combinations. 

For other statistical fluctuation and detector effects of different methods compari- 
son, please refer to figure DOT] for positive charge and figure DOl2l for negative charge. 
They are further separated by r] > and rj < regions and UD/LR directions. To 
check for consistency, we plot the ratios of the {A 2 stat+de1 ) between the positive and 



negative charges and the UD and LR directions in figure 3.13(d) , data points are ob- 
tained by adding- n method. Since the events are divided into hemispheres, the total 
multiplicity of the charged particles are same in UD and LR directions. In principle, 

( A ud, stat+det) and ( A LR, s tat+det) are the same > which is shown with the triangle sym- 
bols for positive charge in 77 > and 77 < regions. The ratios are consistent with 
unity for all centralities. However, the ratios between positive and negative charges, 
shown in circles, are consistent with 98% instead of unity for all centralities. This is 
because the total multiplicity of the positively and negatively charged particle average 
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multiplicities are different due to isospin asymmetry. We collide gold nuclei with net 
charge of positive 79 from each gold ion. In the TPC covered — 1 < 77 < 1 region, the 
total multiplicity of recorded positively charged particles is then roughly 2% larger 
than the multiplicity of negative charges. As we mentioned, the statistical fluctuation 
is approximately proportional to the inverse of the total multiplicity. Then, the ratio 
of between (A 2 +LR>stat+det ) and (A 2 _ LRstat+det ) is off unity by 2%. We also show the 
ratio of positive and negative charges using "1/7V" approximation in black line, which 
is consistent with the adding- 7r method. 

For simplicity and computing efficiency reason, we use adding-7r method for our 
statistical fluctuation and detector effects measurement in this analysis. We have 



calculated the adding- ir results in figure 3.4(c) to demonstrate the stat+det for co- 
variances are consistent with zero. 



3.8 Consistency Check 



In the final result, the charge asymmetry correlations are averaged over different 
pseudo-rapidity regions for variances and covariances, and also for different charge 
combinations for variances. In this section, we check all the individual asymmetry 
correlations to make sure they are consistent with each other. 

Figure 13.141 shows the asymmetry variances with their statistical fluctuation and 



detector effects superimposed in lines. We show {A 2 +UD ) in figure 3.14(a)), {A 2 + LR ) 
in figure 3.14(b)), (A 2 _ UD ) in figure 3.14(c) and (A 2 _ LR ) in figure 3.14(d) with p T 



range of 0.15 < px < 2.0 GeV/c as a function of N part . The red data points and lines 
denote the variances and their statistical fluctuation and detector effects from 7] < 
region, and blue ones are from rj > region. The asymmetry variances are scaled by 
Npart to better show the differences. From the plots we can see all the asymmetry 
correlations are with similar magnitude and centrality dependence. They are largely 
dominated by the statistical fluctuations. All variances in rj < region (red squares) 
are slightly larger than the corresponding variances in 77 > region (blue circles). 
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Figure 3.14. Consistency check of the variances as a function of centrality 
for panel (a): (A^ UD ),(b): (A 2 +jLR ), (c): (A 2 _^ UD ), (d): (A\ LR ) and their 
statistical fluctuation and detector effects. And the dynamic variances 
panel (e): 6(A 2 UD ) and (f): 5(A 2 LR ). The data are RUN IV 200 GeV 
Au+Au collisions with particle pt range of 0.15 < pr < 2 GeV/c. All 
correlations are scaled by N part . 
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This is because of the detector efficiency deficit in r] < region. The estimate of 
statistical fluctuation plus detector effects using adding-7r method is also shown in 
lines in each plot. The red line shows the t] < region, which is also slightly larger 
than the blue line (r] > region) because of the same reason. 

After subtracting the statistical fluctuation plus detector effects, we obtain the 



dynamical variances in UD direction, which are shown in figure 3.14(e), and in LR 



directions which are shown in figure 3.14(f) The dynamical variances show good con- 
sistency between the different pseudo-rapidity regions and different charges. There- 
fore, we take the average over the different combinations to increase the statistics. 
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Figure 3.15. Consistency check of the covariances and their statistical 
fluctuation and detector effect scaled by N part as a function of N part . The 
data are RUN IV 200 GeV Au+Au collisions with particle pr range of 
0.15 <p T < 2.0 GeV/c. 



The covariance consistency check is shown in figure 13.151 The statistical fluctua- 
tion plus detector effects are consistent with zero for all centralities and combinations, 
so we do not have to subtract them. Again we take the average of the different com- 
binations. 
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3.9 Systematic Uncertainties 

In order to study systematic effects, we have studied various conditions and cuts 
made in the analysis. 

3.9.1 Sanity checks 

To show the charge asymmetry correlations are indeed caused by the event-plane 
effect, we use random event-plane to check the asymmetry correlations. The anal- 
ysis procedure is the same, except we generate a uniform distributed number <Pep 
in range of < <Pep < 2tt as the event-plane azimuthal angle. The asymmetries 
and their correlations are then calculated with respect to the randomly generated 
event-plane direction. The results are shown in figure 13.161 where the asymmetry 
correlations are scaled by N part to better show the scale. Because the event-plane is 
random, the asymmetries of UD and LR cannot be distinguished. Thus the asymme- 



try correlations are then identical for UD and LR, which is shown in figure 3.16(a 



We can see that both the same-sign and opposite-sign asymmetry correlations are con- 



sistent between UD and LR. Figure 3.16(b) shows the difference between UD and 



LR correlations. They are consistent with zero for the same-sign and opposite-sign 
correlations, i.e. (A + A^)ud = (A + A_)lr and (Afj D ) = (A 2 LR ). These results give us 
confidence that the charge asymmetry correlations we present are indeed caused by 
the event-plane. The asymmetry correlation differences between the U D and LR are 
indeed correlated to the event initial geometry configuration. 

We also rotate the reconstructed event-plane by ir/4 counterclockwise to study the 
systematic uncertainty. The asymmetries and their correlations are then calculated 
with respect to the rotated event-plane. The results are shown in figure 13.171 After 
rotation, the UD and LR are shifted, so that the new UD and LR hemispheres are 
equally mixed by the original UD and LR hemispheres. As expected, the asymmetry 
correlations are destroyed, UD and LR correlations cannot be distinguished as seen 
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(a) Dynamical asymmetry correlations 



(b) UD — LR asymmetry correlations 



Figure 3.16. Panel (a): Centrality dependences of the charge multiplicity 
dynamical correlations 5{A 2 ) and the opposite-sign correlations (A + AJ) 
for RUN IV 200 GeV Au+Au collisions with respect to random EP. Panel 
(b): The UD — LR correlations of the same-sign and opposite-sign cor- 
relations. The asymmetry correlations are scaled by the number of par- 
ticipants N par t to better show the magnitude. The particle pr range is 
integrated over 0.15 < pr < 2.0 GeV/c for the asymmetry calculation. 
Error bars are statistical errors only. 



in figure 3.17(a) The differences between UD and LR are consistent with zero, and 



show no centrality dependence in figure 3.17(b) 



To estimate the directed flow effects, we separate the event with pseudo-rapidity 
range of \r]\ < 0.5 and 0.5 < |^| < 1.0. Since the directed flow is an odd function 
of rj, and close to zero in \rj\ < 0.5 region, we calculate the asymmetries and their 
correlations from sub-event in \rj\ < 0.5 with respect to the event-plane reconstructed 



from 0.5 < \rj\ < 1.0 region. The results are shown in figure 3.18(a), which are 
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(a) Dynamical asymmetry correlations 



(b) UD — LR asymmetry correlations 



Figure 3.17. Panel (a): Centrality dependences of the charge multiplicity 
dynamical correlations 5{A 2 ) and the opposite-sign correlations (A + A_) 
for RUN IV 200 GeV Au+Au collisions with respect to EP rotated by 
7r/4. Panel (b): The UD — LR correlations of the same-sign and opposite- 
sign correlations. The asymmetry correlations are scaled by the number 
of participants N part to better show the magnitude. The particle pt range 
is integrated over 0.15 < pt < 2.0 GeV/c for the asymmetry calculation. 
Error bars are statistical errors only. 



consistent with using rj < and t] > separated sub-event results in figure 14.21 For 



better comparison, we show the difference between the two in figure 3.18(b) , and also 



the differences between the UD — LR correlations in 3.18(c) The difference is close 
to zero, and does not depend on centrality. 
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Figure 3.18. Panel (a): Asymmetry correlations scaled by number of 
participants N part as a function of centrality for particles from \q\ < 0.5 
region and event-plane reconstructed from particles from 0.5 < \q\ < 1.0 
region. Panel (b): The differences between the correlations from \rj\ < 0.5 
regions and the half TPC sub-event correlations (figure Panel (c): 

The differences between the UD — LR correlations from \rj\ < 0.5 regions 
and the half TPC sub-event correlations. 



3.9.2 Analysis Cut Variations 



The systematic uncertainties can be estimated by varying various STAR standard 
cuts. Firstly, we check the magnetic field polarity uncertainties. We separate the 
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data according to the magnetic field direction, full magnetic field (FF) and reversed 
full magnetic field (RFF). We use C collectively denotes the final correlation values 
of either variances or covariance with any charge or UD / LR directions combinations 
with standard cuts shown in table 13.21 For different magnetic fields, we have the 
correlations Cff and Crff separately. The systematic uncertainty Sbf is then esti- 
mated as |C — Cff\ or \C — Crff\ whichever is larger, where C is the magnetic field 
integrated value. The systematic uncertainties are shown in figure [37191 as red lines. 

Secondly, we check the vertex position v z cut. The standard cut we used in this 
analysis is \v z \ < 30 cm. We then vary the cut to \v z \ < 15 cm to study the difference 
|C|d,|<i5 — C\ which is denoted as the systematic uncertainty S Vx . The results are 
shown in figure 13.191 as blue lines. 

Thirdly, we study the standard DCA cut of DC A < 2 cm. We vary it to DC A < 
1 cm and DCA < 3 cm to study the changes of the correlations. The difference 
\C — Cdca<i\ or \C — Cdca<z\ whichever is larger is then treated as the systematics 
of the DCA cut, and noted as Sdca- The results are shown in figure 13.191 as black 
lines. 

Fourthly, the number of fit points can also contribute to the systematics. We 
vary the cut to nfit < 15 from the standard nfit < 20 and treat the difference 
\C — C n fi t< ic\ as another systematic uncertainty S n fu- The results are shown in 
figure 13.191 as green lines. 

Lastly, the asymmetries from east- and west-side of the TPC may contribute to 
the systematics. We use the asymmetry correlations from east- and west-side of the 
TPC comparing to the average values to estimate the uncertainties. The systematics 
Sew are then defined as \C — CEast\ or \C — Cwest\ whichever is larger. The results 
are shown in figure 13.191 as cyan lines. 

The total systematic uncertainties are obtained by combining all five sources 

S = yj S 2 BF + Sl z + S 2 DCA + Sl fit + S 2 EW . (3.14) 
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Figure 3.19. Systematic uncertainties from difference sources as a func- 
tion centrality in dynamical variances panel (a), variances panel (b) and 
the differences between UD and LR panel (c). Refer the text for details. 



The systematic uncertainty values are detailed in table 13.41 for covariances, table 13.51 
for dynamical variances and table 13.61 for the UD — LR differences. And they are 
shown in figure 14.21 and figure 14.31 with the shaded areas. 



Table 3.4 

Systematic uncertainties of the opposite-sign correlations (A+A-). All the numbers are scaled by the correspond- 
ing number of participants N part . 
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Table 3.5 

Systematic uncertainties of the same-sign dynamical correlations S(A 2 ). All the numbers are scaled by the 
corresponding number of participants N part . 
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Table 3.6 

Systematic uncertainties of the UD — LR correlation differences. All the numbers are scaled by the corresponding 
number of participants N part . 
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4. RESULTS AND DISCISSIONS 

In this chapter, we report the charge multiplicity asymmetry correlation results for 
Au+Au collisions at center of mass energy of ^snn = 200 GeV. We also report 
the asymmetry correlations as function of transverse momentum (pr), event-by-event 
anisotropies (f^ 5 ), wedge size and wedge axis location. We discuss the implication of 
the results and compare with previously published results. 



4.1 Charge Asymmetry Correlations 

It is obvious that the single charge multiplicity asymmetries ((A +> ud), (^—,ud), 
(A +i lr), (A-^lr)) are, by definition, zero because the directions of up (left) and down 
(right) hemispheres are all random from event to event. The up (left) and down (right) 
directions only have relative meaning, and cannot be measured experimentally. The 



data indeed show that they are centered at zero within statistical errors in figure 14.11 
where data are from RUN IV 200 GeV Au+Au collisions. The asymmetries are calcu- 
lated from one side of the TPC tracks with respect to the event-plane reconstructed 
from the other side of the TPC tracks with pt ranges of0.15<pr<2.0 GeV/c for 
both the asymmetry calculation and event-plane reconstruction. The average values 



of the asymmetries as a function of the centrality are shown in figure 4.1(a) for UD 



direction and figure 4.1(b) for LR direction. They are multiplied by the number of 
participants N part . 

We report the variances and covariances of the charge multiplicity asymmetries 
shown in figure 14.21 which are obtained from the Au+Au and d+Au collisions with 
the center of mass energy of a/J/vjv = 200 GeV. The Au+Au data are from RUN 
IV with centrality cuts and N part shown in table 13.11 d+Au data are from RUN 
III, and Glauber model shows the average number of participants for d+Au collision 
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Npart '"'part 

(a) (A) UD (b) (A) LR 

Figure 4.1. Mean value of the charge multiplicity asymmetries of UD di- 
rection in panel (a) and LR direction in panel (b). The data used are from 
Au+Au 200 GeV RUN IV with p T ranges of 0.15 < p T < 2.9 GeV/c for 
both the asymmetry calculations and event-plane reconstruction. Error 
bars are statistical errors only. 

at y / J/v r /v = 200 GeV is 8.31 j28j[3l]. All the charge asymmetry correlations are 
scaled by the number of the participants N part . The dynamics and properties for 
d+Au collisions are very similar to the very peripheral Au+Au and p+p (proton 
to proton) collisions, so we present the d+Au result as a reference. It is believed 
that the QGP can hardly be formed in p+p, d+Au and very peripheral Au+Au 
collisions, hence chiral symmetry cannot be restored in such collisions. Therefore the 
correlations in d+Au and the most peripheral Au+Au collisions are the background 
that is not due to CME/LPV. We present d+Au data to show the trend in the 
correlations is consistent when moving from Au+Au peripheral collisions to d+Au 
collisions. Note that in d+Au collisions, the reaction-plane is not defined because 
there is no anisotropy flow in such small multiplicity collisions. If we apply the same 
algorithm of the second order event-plane reconstruction as shown in section l3.5.1[ the 
event-plane reconstructed is mostly dominated by fluctuation and/or di-jets, similar 
to the peripheral Au+Au collisions. 
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Figure 4.2. Centrality dependence of the charge multiplicity dynamical 
correlations S(A 2 ) and the opposite-sign correlations (A+A_) for RUN 
IV 200 GeV Au+Au collisions. The asymmetries are calculated between 
hemispheres from a half event divided by event-plane {UD) reconstructed 
from the other half event, and the plane perpendicular to the event-plane 
(LR). The asymmetry correlations are scaled by the number of partici- 
pants Npart to better show the magnitude. The particle pt range is inte- 
grated over 0.15 < pt < 2.0 GeV/c for both the asymmetry calculation 
and event-plane reconstruction. The upper and lower shaded areas illus- 
trate the half size of the systematic uncertainties on (A + AJ) and (A 2 ) 
respectively in UD or LR direction whichever is larger of each centrality. 
The error bars are smaller than the symbols. 



The dynamical asymmetry variances (same-sign correlations) are shown in fig- 
ure 14.21 as square symbols. The solid squares are the dynamical variance in UD 
direction (5{A 2 JD )), and open squares are the dynamical variance in LR direction 
(S(Aj jR )). A positive variance 5 (A 2 ) indicates broadening of the single asymmetry 
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distributions of A + and A_ due to any dynamical processes beyond statistical fluctu- 
ation, such that the same-sign particles are emitted more likely in the same direction. 
On the other hand, a negative variance indicates narrowing of such distributions due 
to symmetric (back-to-back) correlations. From figure 14.21 in peripheral collisions, 
both 5(Af rD ) and 5{A 2 LR ) are positive, which suggests that the same-sign particles 
are more likely emitted in the same direction (i.e., they have small angle correlation) 
within one unit of pseudorapidity, —1 <?7<OorO<?7< 1. When moving toward 
more central collisions, both of the variances become negative. Same-sign particles 
are more preferentially emitted back-to-back, in other words, more symmetrical in 
central collisions than in peripheral. And the small angle correlation is stronger in 
UD direction than that in LR direction, i.e. 5(A 2 JD ) is larger than 5(A\ R ) for all 
centralities. 

Figure 14721 also shows the covariances (A + A_)ud in solid circles and (A + A_)lr 
in open circles. The correlations between the positive and negative asymmetries are 
large and positive for both UD and LR directions. They are on the order of about 
1CT 3 before scaled by N part , suggesting the individual asymmetries are as large as 
a few percent, comparable to the elliptic flow (vz) magnitude. As discussed, posi- 
tive correlations suggest that the positive and negative charged particles tend to be 
emitted in the same direction (small angle correlation). The small angle correla- 
tion is stronger in the UD than that in LR direction in most centralities. And the 
correlations are much stronger in the opposite-sign than the same-sign. 

We also show the same-sign and opposite-sign correlations in d+Au collisions 
in the leftmost symbols. Those correlations follow the peripheral collisions of the 
Au+Au collision, suggesting a smooth trend. 

Figure 14.21 depicts the following picture. d+Au collisions and very peripheral 
Au+Au collisions fall in the same trend. Particles, within one unit of pseudorapid- 
ity (—1 < rj < and < r\ < 1), are preferentially emitted in the same direction 
(with small angle correlation), no matter same-sign pairs or opposite-sign pairs. The 
magnitude of the small angle correlation is, however, stronger in the opposite-sign 
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than that in same-sign pairs, and is stronger in the out-of-plane (UD) than in-plane 
(LR) direction regardless of charge combinations. From mid-central to central col- 
lisions, the same-sign pairs are preferentially back-to-back, while the opposite-sign 
pairs are still preferentially aligned in the same direction and more so than in periph- 
eral collisions. The opposite-sign pair small angle correlation is always stronger in the 
out-of-plane (UD) than in-plane (LR) direction. On the other hand, the back-to-back 
emission tendency of same-sign pairs is weaker in the out-of-plane (UD) that in-plane 
(LR) direction. It is interesting to see that same-sign pair correlations show different 
behavior in peripheral collisions and most central collisions, which implies different 
mechanisms dominate in each situation. 

The systematic uncertainties are shown in the shaded areas for the covariances 
(upper band) of UD or LR direction whichever is larger, and for the variances (lower 
band) of UD and LR whichever is larger. The systematic uncertainties are negligible 
compare to the magnitude of the correlations. 

In order to investigate the possible CME/LPV effect, we study the difference 
between UD and LR asymmetry correlations of the variance (A(A 2 ) = 5(A 2 ID ) — 
5(A 2 LR )) and covariance (A(A + A_) = (A + A_)ud — (A + A_) lr). The charge inde- 
pendent background not related to event-plane cancels by taking the difference. For 
variances, we take the difference between the dynamical correlations, where statistical 
fluctuation and detector effects have been subtracted according to different charges 
and rj regions respectively. 

Figure 14731 shows the centrality dependence of A{A 2 ) and A{A + A_) scaled by the 
number of participants N part . Also shown in curves are the A(A 2 ) and A(A + AJ) 
values with event-plane resolution correction. It is done by linear extrapolation to 
unity of the correlations if higher-order harmonic terms can be ignored. 

The shaded areas are the systematic uncertainties of A(A 2 ) in upper band and 
A(A + A_) in lower band both multiplied by N part . 

From figure H~3l we can conclude the following. First, same-sign variance is positive 
A(A 2 ) > (5(A 2 JD ) > S(A 2 LR )) in all centralities. The UD asymmetry distribution 
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Figure 4.3. The correlation differences between out-of-plane UD and 
in-plane LR of the same-sign A(A 2 ) = (Afj D ) — (A 2 LR ) and opposite-sign 
A(A + A_) = (A + A-)ud — (A + A_) lr. Both correlations are multiplied by 
Npart- The particle pr range of 0.15 < pr < 2.0 GeV/c is used for both 
asymmetry calculation and event-plane reconstruction. The leftmost data 
points are from d+Au data. The curves are the linear-extrapolated val- 
ues of the A(A 2 ) (blue) and A(A + A_) (red) corresponding to the perfect 
event-plane resolution of unity. The shaded bands are the systematic un- 
certainties of A (A + A_) (upper band) and A (A 2 ) (lower band). Statistical 
errors are smaller than the symbols. 



is always broader than the LR distribution. The dynamical same-sign correlations 
S{A 2 ) are the width of single asymmetry correlation beyond the statistical fluctuation. 
More small angle pairs will increase the dynamical correlation, while more back-to- 
back pairs will decrease it. Thus, there are more small angle pairs in the out-of-plane 
(UD) direction than that in the in-plane (LR) direction, or, more back-to-back pairs 
in-plane than out-of-plane. Or, it is possible both cases are true. 
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Second, opposite-sign covariance is positive A (A + A_) > ((A + A-)jjd > (A+A-)lr) 
in all centralities except the most peripheral bins. Similar to the discussion above, 
more small angle opposite-sign pairs are emitted in out-of-plane (UD) direction than 
that in in-plane (LR) direction. 

Third, the magnitude of opposite-sign correlation A(A + A_) is small relative to the 
correlations ((A + A_)ud and (A + A_) lr) themselves. This indicates that the majority 
of the strong background correlations between opposite-sign pairs are unrelated to the 
reaction-plane. The UD — LR difference, which is related to the reaction plane, is a 
small effect compared to the background. 

Lastly, A(Afj D ) and A(A + A^) have similar centrality dependence. They both 
increase with centrality from most peripheral collisions, and reach a maximum in 
medium central collisions. Then, they decrease with centrality to almost zero in 
most central collisions. While the opposite-sign A(A + A_) starts slightly negative 
in very peripheral collisions, which is also true for d+Au collisions for both same- 
sign and opposite-sign correlations. It is mainly due to non-flow effect, such as di- 
jets, dominating the EP reconstruction. The decrease towards most central collisions 
is most likely due to the lower EP resolution and the dilution effect by the large 
multiplicities in the most central collisions. In most central collisions, the low EP 
resolution cannot distinguish between up and down, or left and right. Thus the 
difference between U D and LR vanishes. 

CME/LPV predicts that charge separation is along the system angular momentum 
direction [7J[9j[l0]. In the absence of the medium effect, the charge asymmetry cor- 
relations between light quarks across the event-plane will carry to the final hadronic 
state. Thus the final measured charged particle correlations are expected to have 
the charge separation effect [S|. CME/LPV expects that more negatively charged 
particles going down (up) across the reaction-plane, if more positively charged par- 
ticles going up (down) across the reaction-plane. In other words, positive and neg- 
ative charges are anti-correlated in UD direction. On the other hand, CME/LPV 
effect does not contribute to the LR direction. Therefore the covariance in UD di- 
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rection should be smaller than that in LR direction, (A+A-)ud < (A+A-)lr, or 
A(A + A_) < 0. Since CME/LPV generates additional correlations in UD direction 
rather than LR direction, one should expect the variance in UD is larger than that in 
LR direction. This can be written in the language of charge asymmetry correlations 
as 5(Al tUD ) > 5(Al LR ), or A(A 2 ) > 0. 

CME/LPV effect has to occur in the early time of the collision while the low mass 
quarks are deconfined, chiral symmetry restored, and the magnetic field is large [10] . 
The correlations due to CME/LPV could be modified by the interactions with the 
hot dense medium. Effect of medium interactions has been observed in the jet-like 
two-particle correlations with strong back-to-back correlation suppression at high pt 
and enhancement at low px [35H37] . One may expect similar suppression effect on the 
back-to-back charge asymmetry correlations, which are the opposite-sign correlations 
across the reaction-plane. 

Data shown in figure 14.31 seem to suggest that the same-sign correlation is consis- 
tent with CME/LPV expectation that additional correlations broaden the asymme- 
try distribution in UD direction. However, the dynamical variances of 5{A 2 JD ) and 
5(A 2 LR ) in figure H~2l show that they are negative in mid-central to most central col- 
lisions. In other words, the same-sign particles are more likely emitted back-to-back, 
which is inconsistent with CME/LPV expectation. The opposite-sign correlations 
are also inconsistent with the expectation that the opposite-sign pairs are strongly 
correlated in the same direction, and the effect is stronger in UD than LR. It appears 
that CME/LPV alone cannot explain the data. 

As we discussed in section II. 3[ the magnitude of the opposite-sign correlation 
expected from CME/LPV is negative around ~ 10~ 5 to ~ 10~ 4 , and ~ 10~ 3 to 
~ 10 -2 after multiplied by the number of participants [THTO]. The covariances in UD 
and LR directions are significantly larger than the expected value by 1 to 2 orders and 
they are positive, which suggests that they are largely dominated by the event-plane 
unrelated background. After taking the difference between UD and LR, it has the 
order of 10~ 2 for mid-central collisions but still in positive. 
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It has been suggested that the charge separation effect can be invoked by QCD 
processes without the need of CME/LPV to create such asymmetry correlations [38J. 
However, the effects are estimated to be orders of magnitude smaller than our observed 
charge asymmetry correlations. We note that such estimates have large uncertainties 
due to uncertainty in various effects, such as the effect of medium interaction [T4"t[T5]. 



4.2 Charge Asymmetry px Dependence 

We study the px dependence of the asymmetry correlations which has been shown 



in figure 4.4(a) for 0-20%the most central collisions, figure 4.4(b) for 20-40% mid 



central collisions, and figure 4.4(c) for 40-80% peripheral collisions. The px ranges 
for particles used in asymmetry calculation are 0.15-0.35, 0.35-0.5, 0.5-0.75, 0.75-1.0, 
1.0-1.25, 1.25-1.5, and 1.5-2.0 GeV/c. The event-plane is still reconstructed from 
particles with 0.15 < px < 2.0 GeV/c. All dynamical variances and covariances in 
U D and LR are positive at very low px for all centralities. The dynamical variances 
S(A 2 ) shown as squares drop rapidly to zero with px up to 1 GeV/c (soft particles) 
for all centralities. The 5(Aj jR ) becomes negative at px > 1 GeV/c, and the 5(Afj D ) 
is around zero, slightly negative at px > 1 GeV/c for all centralities. The same-sign 
correlation px dependence results indicate that, small angle pairs dominate at low-pr 
{px < 1 GeV / c) . The tendency for back-to-back emission of same-sign pairs increases 
with increasing px- 

On the other hand, opposite-sign correlations (A + A^) shown as circles remain 
relatively constant with px up to 1 GeV/c, and then rapidly increase with increasing 
Px- The results indicate that, opposite-sign pairs are emitted in the same direction at 
low-pr, and the small angle correlation increases strongly with px > 1 GeV/c for both 
out-of-plane (UD) and in-plane (LR) directions. The correlations are qualitatively 
similar for all centralities. 

Figure 14.51 shows the UD — LR of same-sign and opposite-sign correlations as a 



function of px for the most central in figure 4.5(a), mid-central in figure 4.5(b) and 
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Figure 4.4. The px dependence of the charge asymmetry dynamical 
correlations 5 (A 2 ) and the opposite-sign charge asymmetry correlations 
{A + A^} for RUN IV Au+Au 200 GeV collisions. The most central 0-20% 
collisions are shown in panel (a), medium central 20-40% collisions are 
shown in panel (b), and the peripheral 40-80% collisions in panel (c). The 
asymmetries are calculated between hemispheres separated by event-plane 
(UD) and the plane perpendicular to the event-plane (LR). The particle 
Pt range of 0.15 < pt < 2.0 GeV/c is used for event-plane reconstruction. 
Error bars are statistical only. 



peripheral collisions in figure 4.5(c) Both same-sign and opposite-sign correlations 
grow with increasing px, and show similar p? dependence. The same-sign corre- 
lation is stronger than the opposite-sign correlation in mid-central collisions, while 
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Figure 4.5. The dependence of the charge asymmetry correlation 
differences between UD and LR of same-sign A{A 2 ) = 5(A 2 JD ) — 5(A 2 - yR ), 
and opposite-sign A(A + A_) = {A + A^)uo — (A + A^)lr. The event-plane 
is constructed from particles pt range of 0.15 < pt < 2.0 GeV/c. Error 
bars are statistical. 



the difference between same-sign and opposite-sign correlations is not obvious in the 
most central and peripheral collisions. 

The CME/LPV effect is non-perturbative and subjected to soft particle produc- 
tion correlation. The limit of the transverse momentum p? range for such effect being 
experimentally tested at center of mass energy ^/snn = 200 GeV is estimated smaller 
than 1 GeV/c [TO], which may be subjected to radial flow effect. Thus, CME/LPV 
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expects charge separation at low-p^ and little separation at high-p^. However this is 
not obvious in the data due to limited statistics. 

It is worthwhile to mention that the charge asymmetry correlations for each p? 
range, in this section, are calculated only within the px bin. Meanwhile, the centrality 
dependence of figure H~2l and figure H~3l is calculated from all particles within 0.15 < 
Pt < 2.0 GeV/c. As such, the data shown in the previous section cannot be readily 
obtained from those in figure 14.41 of the corresponding centrality bin. 

4.3 Event-by-Event Anisotropy {v^ s ) Dependence 

The physics mechanism for such event-plane dependent charge asymmetry corre- 
lations is unclear. There are alternative models other than CME/LPV effect, such 
as, general cluster particle correlations with anisotropies that could generate sizeable 
difference between in-plane and out-of-plane particle correlations [39]. Or, the mo- 
mentum conservation and local charge conservation together with elliptic flow could 
yield event-plane dependent correlations that differ between same- and opposite-sign 
pairs [TUHO]. Path-length dependent jet-quenching effect could be another mecha- 
nism [37J|39llJT] . For most of these, one qualitatively expects event anisotropy depen- 
dence. 

The similar trends of same-sign and opposite-sign charge asymmetry correlations 
in figure H~3l are qualitatively consistent with that of elliptic flow as a function of cen- 
trality [18]. The asymmetry correlations have maxima at mid-central collisions, and 
drop in central and peripheral collisions for both same-sign and opposite-sign corre- 
lations. Motivated by these considerations, we examine the event shape dependence 
of the asymmetry correlations. 

We investigate the dynamical variances 5(A 2 ), covariances (A+A-), and their 
difference between out-of-plane (UD) and in-plane {LB) A(A 2 ) and A(A + A_) as a 
function of event-by-event azimuthal anisotropy {v^" 3 ) of high-py {px > 2 GeV/c) and 
\ow-pr {pr < 2 GeV/c) particles respectively. The high-p^ anisotropy dependence 
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may be sensitive to the jet-quenching effect where the high-pr particles are suppressed 
in out-of-plane direction [HT|[57] . The \ow-px anisotropy dependence may characterize 
the bulk event shape of the charge combinations. 

The event anisotropy is defined as v 2 bs = (cos 2(0 — i/jep)) with low-p^ and high- 
Pt particles. For low-pr anisotropy (cos 2(0 — ipEp)), the average is taken from all 
particles used in asymmetry calculation (within one half of the TPC), and the event- 
plane is reconstructed from the other half of the TPC. The v^ t<2 GcV / c denotes the 
low-pr anisotropy. For high-pr anisotropy, in order to increase statistics, particles 
with pt > 2 GeV/c from the entire event (|?7| < 1) are used for the anisotropy 
calculation, while the EP is reconstructed from different sides of the TPC. We use 
u 2pt>2 GeV/c t° stand for high-pr anisotropy. 

A positive v 2 bs indicates an event with more particles of the interest emitted in- 
plane, i.e. the event is elongated in the event-plane direction. On the other hand, a 
negative v 2 bs indicates more interest particles emitted out-of-plane, i.e. the event is 
elongated perpendicular to the event-plane. A zero v 2 bs means the event is spherical, 
i.e. the emitted particles are isotropic in azimuth. We use term "spherical" to refer 
events with such shape. 

Figure H~6l shows the asymmetry correlation results in mid-central 20-40% central- 
ity Au+Au 200 GeV RUN IV collisions as a function of event-by-event anisotropy 
for high-p T v°*p T>2 GcV/c . The high-p T anisotropy {v 2 p T>2 GeV/c ) dependence of same- 



sign dynamical variances S(A ) is shown in figure 4.6(a), and the opposite-sign co- 



variances (A + AJ) is shown in figure 4.6(b). No significant high-p^ event-by-event 
anisotropy dependence is observed of same- and opposite-sign correlations in both in- 
plane and out-of-plane directions, which may suggest that the path-length dependence 
jet-quenching has little effect on the charge asymmetry correlations. The differences 



between the UD and LR of the variance and covariance are shown in figure 4.6(c) 
Both A(A 2 ) and A(A + A_) are positive and flat over a large range of the high-pr 
anisotropy. This is consistent with the finding of medium modified jets and initial 
state fluctuations in 
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Figure 4.6. High-p^ event-by-event anisotropy v^ s dependence of (a) 
the dynamical charge asymmetry variances, (b) asymmetry covariances, 
and (c) their differences between UD and LR hemispheres, in mid-central 
20-40% Au+Au 200 GeV collisions from RUN IV. The asymmetries are 
calculated from all particles from one side of the TPC tracks with respect 
to the event-plane reconstructed from the other side of the TPC tracks. 
Both the particles used in asymmetry calculation and event-plane recon- 
struction are with range of 0.15 < px < 2.0 GeV/c. The asymmetries 
are correlated to the high-pr event-by-event anisotropy , which is cal- 
culated from particles with pt > 2.0 GeV/c in the entire TPC to increase 
statistics. Error bars are statistical only. 



Figure 14.71 shows the low-pp anisotropy dependence of the charge asymmetry cor- 
relations. Different from high-pr anisotropy, significant v^ 8 dependence is observed 



in the dynamical variances 5(A ) in figure 4.7(a) The out-of-plane variance 5(Afj D ) 
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increases with low-pr v^ s , while the in-plane variance 5(A 2 LR ) decreases with v^ s . 
The different trend results in a strong vf" 3 dependence of the difference between U D 



and LR (A{A )), which is shown in figure 4.7(c) in squares. Some low-p^ v% de- 
pendence is observed in the covariances (A + A— ), but is significantly weaker than the 
variances. However, we can still see opposite trend in the covariances, that a slightly 
decreasing trend for low-p^ mUD and increasing trend in LR as shown in figure 



4.7(b) The difference of the covariance A(A + AJ) is shown in figure 4.7(c) in circles, 
which decreases with increasing v^ s . It is interesting to see that the same-sign and 
opposite-sign UD — LR correlations have different trend with low-pr v^ s . 

It is important to point out in figure 4.7(c) that, A(A 2 ) and A(A + A^) cross at 
v^ s ~ 0, and the crossing value is positive. For events with v 1 ^" 3 ~ 0, no significant 
charge difference is observed between same-sign and opposite-sign pair correlations, 
at which the particle distribution is isotropic. We know that the average v^ s over the 
entire event sample is positive due to the elliptic flow, table 13.11 The different trend 
of the variance and covariance as a function of vf" 3 leads to the integrated values of 
variance and covariance to diverge. In other words, after integration, the UD — LR 
variance is larger than the covariance, i.e. A{A 2 ) > A{A + A_) for all centralities at 

( V 2,pr<2GeV/c) > 0- 

There is centrality dependence of the magnitudes of the asymmetry correlations 
as a function of v% bs . We show the central 0-20% Au+Au 200 GeV v^ s asymmetry 
correlation dependence in figure IA. 131 and the peripheral 40-80% collisions in figure 
IA.14l in the appendix. The features are qualitatively similar as the mid-central colli- 
sions. For high-pr v% bs , the dynamical variances and covariances are consistent with 
constant for all centralities. The difference between UD and LR does not depend 
on high-py v% bs . However, the dynamical variances and covariances strongly depend 
on 1ow-pt v% bs , and the difference between UD and LR shows opposite trend of the 
variance and covariance in all centralities. The magnitudes are different though. 

To further remove possible short range correlation between the particles used in 
asymmetry calculation and the event-plane reconstruction, we choose particles with a 
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Figure 4.7. Low-p^ event-by-event anisotropy v^ s dependence of (a) 
the dynamical charge asymmetry variances, (b) asymmetry covariances, 
and (c) their differences between UD and LR hemispheres, in mid-central 
20-40% Au+Au 200 GeV collisions from RUN IV. The asymmetries are 
calculated from all particles from one side of the TPC tracks with respect 
to the event-plane reconstructed from the other side of the TPC tracks. 
Both the particles used in asymmetry calculation and event-plane recon- 
struction are with pt range of 0.15 < pr < 2.0 GeV/c. The asymmetries 
are correlated to the low-pr event-by-event anisotropy v^ s , which is cal- 
culated from the same particles used for the asymmetries. Error bars are 
statistical only. 



large pseudo-rapidity gap between the two particle sets for event-plane reconstruction 
and the charge multiplicity asymmetry correlation. Instead of separating event by 
— 1 < rj < and < r\ < 1, we separate the event with —1.0 < r\ < —0.5 and 
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Figure 4.8. Event-by-event dependence of the dynamical variances 
(top panel), covariances (middle panel) and UD — LR (bottom panel) of 
RUN IV 20-40% 200 GeV Au+Au data. Particles used for asymmetries 
and EP reconstruction are separated by one unit in 77. Asymmetries are 
correlated to 1ow-pt (left column) and high-p^ (right column) ujj . Error 
bars are statistical only. 
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0.5 < 7] < 1.0, such that the two sub events have a one unit pseudo-rapidity gap, 
which is the usual soft particle bulk correlation span in pseudo-rapidity space. 

We show similar v^" 1 dependence results of mid-central Au+Au 200 GeV data in 
figure H~8l Low-p<r v^ s (jp? < 2.0 GeV/c) is calculated from the same particles used 
for the asymmetries. High-p^ is calculated from all particles with p? > 2.0 GeV/c 
to increase statistics. The corresponding UD — LR correlation with rj gap is shown 
in the bottom panel in figure 14.81 The result is consistent with what we observed 
without rj gap but with larger error bars because of the limited statistics. The central 
and peripheral collision results are shown in figure IA. 151 and figure IA.16I They are 
also qualitatively consistent with the mid-central result. 

To even further remove the short range correlation, we use the first order event- 
plane reconstructed with the ZDC-SMD detector, which is \rj\ > 6 away from the 
main TPC in pseudo-rapidity direction. The large distance in pseudo-rapidity space 
removes short range correlation almost completely. We calculate the charge asymme- 
try correlations and the event-by-event v^ s with respect to the first order ZDC-SMD 
event-plane using RUN VII Au+Au 200 GeV data. We show the result in figure 
14.91 which has very good consistency with the second order event-plane result. The 
central and peripheral asymmetry correlations vs v^ s results are shown in figure lA~T7l 
and figure IA.18[ where similar dependence is observed. 

One other suggestion is to evaluate charge separation effect in the top 2% most 
central collisions [19]. In such events, the magnetic field generated by the wounded 
nuclei is moderate with less fluctuation over time [12]. The fluctuation of event-by- 
event vf s is also small in the very most central collisions. Thus CME/LPV should 
be relatively constant, and not depend on the event shape v^ s . To select such top 
2% central events, we use the ZDC triggered RUN IV 200 GeV Au+Au events, which 
triggers on the top 12% most central collisions. In addition, we require the sum of the 
ADC signal from ZDC-SMD detectors to be less than 78, which will roughly select 
the top 2% top central collisions. We apply the standard STAR quality cuts as we 
do in the minimum-bias triggered data, and the corresponding acceptance correction. 
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Figure 4.9. Asymmetry correlations vs event-by-event v^ s of RUN VII 
200 GeV Au+Au mid-central 20-40% collisions. The particles used for 
asymmetry calculation are from half side of the TPC with respect to 
the first order event-plane reconstructed from ZDC-SMD. Error bars are 
statistical only 



There are in total 5.5 million events of the top 2% central collisions in the final data 
sample. Our results are shown on figure I4.10I Similar features of the event-by-event 
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Figure 4.10. Asymmetry correlations vs event-by-event v^ s of RUN 
IV 200 GeV Au+Au ZDC triggered top 2% most central collisions. The 
particles used for asymmetry calculation are from half side of the TPC 
with respect to the event-plane reconstructed from another side of the 
TPC. Error bars are statistical only. 



anisotropy dependence are seen in the top 2% most central collisions as for 0-20% 
centrality (figure lA.13j) . but the error bars are large due to limited EP resolution. 
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To summarize, in this section we presented the asymmetry correlations in Au+Au 
200 GeV 20-40% mid-central centrality as a function of event-by-event v^ s defined 
using the second order event-plane with and without r] gap, first order event-plane 
from ZDC-SMD detector, and we also checked the results against those for the top 
2% most central collisions. All the results are consistent with each other, which gives 
us confidence that the v^ s dependence is indeed the main cause of the observed final 
state charge asymmetry correlations. 



4.4 Wedge Size and Location Dependence 

All the results we have presented above are charge multiplicity asymmetries ob- 
tained from hemispheres divided by event-plane (UD) and the plane perpendicular 
to the event-plane {LB). The results show the asymmetry correlations, either UD 
or Li?, which are calculated from the same set of particles but divided into different 
sub groups. In this way, statistical correlation and detector effect are mostly canceled 
between UD and LB. However, using the same set of particles will not be sensitive 
enough to distinguish the fine angular charge asymmetry structure. For example, the 
hemisphere study cannot tell us whether the charge separation happens in-plane or 
out-of-plane, and how the charge asymmetry correlations evolve from in-plane {LB) 
to out-of-plane (UD) direction. As we discussed in previous chapter, one of the ad- 
vantages of these observables is that we can vary the azimuthal opening angle A0 W 
in which the analysis particles are included. By restricting the azimuthal range, we 
can study the charge separation as a function of the wedge size and location. The 
wedge size and location study will give us further insight into the fine structure of 
the charge separation. 

We present the wedge size dependence of the mid-central 20-40% 200 GeV Au+Au 
collision data in figure STD for A0 W = 15°, 30°, 45°, 90°, where A0 W = 90° is identical 



to the hemisphere analysis. Figure 4.11(a) shows the charge dynamical asymmetry 
variances S (A^ i±A<h ) , 5(Ag :±A(f) J, and covariances (A + A_) 90 o± A( j >w , (A+A_) o °±A0w 
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(b) UD — LR correlations vs wedge size 



Figure 4.11. The wedge size dependence of charge asymmetry correla- 
tions in panel (a) and their differences between out-of-plane (UD) and 
in-plane (LR) correlations A(A\ l i) and A(A + A_) a</> w in panel (b). Data 
are from 20-40% centrality RUN IV 200 GeV Au+Au collisions. The par- 
ticle pt range of 0.15 < pr < 2.0 GeV/c is used for both asymmetry 
calculation and event-plane reconstruction. Error bars are statistical. 



as a function of the wedge size A<fr w . Note that the statistical fluctuation and detector 
effects are subtracted for each charge and wedge separately in the variances. Both 
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covariance correlations increase with decreasing wedge size A0 W . The dynamical 
variances seem to show similar increasing trend with decreasing A0 W . The smaller 
the wedge size, the stronger the correlations are. This could possibly suggest that the 
charge asymmetry correlations are local, and the major effect is in the opposite-sign 
correlations. 



Figure 4.11(b) shows the difference of the asymmetry correlations between out-of- 



plane and in-plane directions, for variance 

A(A 2 M J = 5(Al QO±M J - 5(Al 0±A<p J (4.1) 

and covariance 

A(A + A^) A ^ = (A + A_) 90 o± 5 ^ - (Ah4_)o°±**w. (4.2) 

Both of them increase with decreasing wedges size, and they show similar trend. The 
difference between these two seems to disappear when A0 W decreases. 

So far, we show the charge asymmetry correlations focusing only on in-plane 
and out-of-plane directions and the difference between them. It is important to 
understand how the correlations evolve from in-plane to out-of-plane. To do so, 
we study the correlation for back-to-back wedges with fixed wedge opening angle 
of 30° (A0 W = 15°), and vary the wedge axis location from in-plane (0 W = 0°) to 
out-of-plane (0 W = 90°). Figure 14.121 shows the dynamical variance and covariance 
of the 30°-size back-to-back wedges in different locations relative to the event-plane 
(0 W = 0°, 30°, 60°, 90°). The data are from mid-central 20-40% centrality 200 GeV 
Au+Au collisions, same as figure 14.111 Both the dynamical variance and covariance 



increase from in-plane to out-of-plane, which is expected from figure 4.11(a), and 
both are modulated by a cosine function (~ cos(20 w )) as shown in the dashed lines. 

We also show the wedge size and location study of the central collisions in figure 
IA. 191 and peripheral collisions in figure IA~.20I Although the magnitudes are different, 
they are qualitatively similar. 
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Figure 4.12. Charge multiplicity correlations as a function of the wedge 
location W of 20-40% centrality 200 GeV Au+Au collisions. The wedge 
opening angle is 30° (A0 W = 15°), and the back-to-back wedges are located 
at 0°, 30°, 60° and 90° relatively to the event-plane. The curves are the 
cosine modulation a x (1 + 26cos(20 w )) to guide the eye. The particle 
Pt range of 0.15 < px < 2.0 GeV/c is used for both the asymmetry 
calculation and event-plane reconstruction. 



4.5 Discussion 

In this section, we first make connections between our observables and the three- 
particle correlators used in previous STAR analysis. Then we discuss the implication 
of our charge asymmetry correlation results. We show the additional information 
from our results of charge separation directions and the event-by-event anisotropy 
analysis. Through the study of event-by-event anisotropy, we may have a further 
understanding of the charge separation signal and/or background. 

4.5.1 Connection to Three-Particle Correlators 

Our study of the charge multiplicity correlations is motivated by the CME/LPV. 
To observe the effect, chiral symmetry restoration and large magnetic field are re- 
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quired [6]. The hot QCD matter is a particle environment to test the effect. CME/LPV 
predicts that light-quark electric charge separates along the direction of the magnetic 
field, which is the same as the system's orbital angular momentum direction. The 
effect yields charge separation in the final state hadrons, such as pions [TJ. 

Previous STAR published result on three-particle correlator (cos (a + (5 — 2c)) 
[iniCnJlEn], we have measured a negative same-sign correlator, and a close to zero, 
even slightly negative opposite-sign correlator. The positive same-sign correlator is 
qualitatively consistent with the local parity violation expectation. However, the close 
to zero opposite-sign correlator seems inconsistent with the naive expectation from 
CME/LPV. Whether the charge separation can survive the hydrodynamic evolution 
to the final state is still an open question. A recent study shows that a large percentage 
of charge separation in the initial state is needed for charge separation to be detected 
in the final state [2] , which suggests the in- medium interaction can strongly modify 
the charge separation correlations. As a result, the back-to-back correlated pairs are 
quenched, because at least one of the pair quarks would be affected by the in-medium 
interactions. Such in-medium effects could qualitatively explain the close to zero 
opposite-sign correlator observed at STAR pTj |[T6l[20] . 

STAR has reported the first measurement of three-particle azimuthal correlator 
to search for CME/LPV in relativistic heavy ion collisions [16l|20]. The correlators 
were introduced in [TS], and can be measured as 



where the average in the last step is taken over event sample, and a and f3 repre- 
sent positive or negative charges. t>2 iC is the elliptic flow of particle c which serves to 
measure the reaction plane, ipnp stands for the reaction plane angle, 5 a p is the Kro- 
necker delta function. In equation I4.3[ it is assumed that the event-plane dependent 
background is negligible, which is not necessarily true. 



(COS(0 Q + (p/3 - 2tp RP )) « (COS(0 Q + - 2<f) c ))/v 2> 




(4.3) 
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Our UD-LR asymmetry correlations A(A|) = 5(A| jl7 _ D )-(J(A| )L _ R ) and A(A + A_) 
(A + AJ) UD — (A + A-)lr are related to the above three-particle azimuthal correlators, 
but with significant differences. The charge multiplicity asymmetries can be expanded 
into Fourier series of step-function in <p — ipEP (particle azimuthal angle relative to 
the event-plane) as 

4 ^ ^ sin[(2n + l)((p±,i - *Pep)} 

A± > UD -^2^2^ 2^Ti ' 

± i=l n=0 

4 ^ ^ cos[(2n + l)(0 ±|i - tp E p)] , AA s 

A ±^ = ^2^2^ ^\ • ( 4 - 4 ) 

= i=l n=0 

The difference between U D and LR is then 

4; 2 

7T 



1 ^v-^' 9 cos[(2n + l){(f> a ,i - iPep) + (2m + - ^ep)] 



^ (2n + l)(2m+ 1) /' 1 ' ' 

a p i,j=0 n,m=0 \ < J\ * J J 

where a and /3 represent positive or negative charges. The asymmetry correlation 
differences contain all possible harmonic terms, including cross terms. While the 
three-particle correlators in equation 14.31 contain only the first order terms. The 
correlators are measured in terms of azimuthal angle with respect to reaction-plane 
angle ipnp, while the asymmetry correlations are measured using azimuthal angle 
relative to event-plane angle ipEP- Thereby, our asymmetry observables are affected 
by event-plane resolution. Hence, the asymmetry correlation observables are related 
to three-particle correlators but are essentially different. 

To gain more insight into the relationship and the differences, we compare the two 
different observables in figure 14.131 In the figure, the charge asymmetry correlation 
LR — UD and the three-particle correlators (cos(0 a + <pp — 2ipnp)) are shown as a 
function of centrality (N part ). The charge asymmetry correlations are shown in solid 
symbols, and the three-particle correlators are shown in open symbols, where "LS" 
stands for like-sign (same-charge), and "US" stands for unlike-sign (opposite-charge). 
The three-particle correlators and LR — UD correlations are scaled by N part . The 
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Figure 4.13. The solid data points are asymmetry correlation LR — UD 
differences of variances A(A 2 ) = 5(A 2 LR ) — 5{A^j D ) (blue) and variances 
A(A + A_) = {A + AJ) LR — (A + A_) UD (red), scaled by the number of partic- 
ipants N part times the scaling constant (7r/4) 2 , and divided by event-plane 
resolution l/e^p . Also shown are the three-particle correlators in open 
symbols, (cos^q,-!-^ — 2(p c ))/v 2tC ~ (cos((p a + (t>i3 — 2ip R p)) of the like-sign 
(blue) and unlike-sign particle pairs (red). The asymmetry correlations 
and the correlator particle a and /3 are calculated with particles from one 
side of the TPC rj < (77 > 0). The event-plane for the asymmetry corre- 
lations and particle c of the three-particle correlators are from the other 
side of the TPC r] > (77 < 0). Particle p T range of 0.15 < p T < 2.0 GeV/c 
is used for asymmetry calculation, correlator calculation and event-plane 
reconstruction. Error bars are statistical only. 



LR — UD correlations are also scaled by a constant (vr/4) 2 according to equation 14.51 
and divided by the event-plane resolution in order to make direct comparison. To 
be consistent, the correlators are calculated using the same set of particles as the 
asymmetry correlations for a and /3, i.e. in the same rj region. Particle c is used from 
the same set of particles to reconstruct the event-plane for asymmetry correlations, 
which is from the other side of the 77 region to avoid self- correlation. The three- 
particle correlator results shown in figure 14.131 are not identical to those in published 
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paper [20], where particle a, (3 and c are from the TPC with pseudo-rapidity range 

Of -1 < 7] < 1. 

From equations 14.31 and I4.4[ we know that the three-particle correlators contain 
only one out of the infinitely many of harmonic terms in the difference of UD and 
LR asymmetry correlations. The same-sign correlation 5(A 2 LR ) — S(A 2 ID ) in figure 
14.131 is comparable to the three-particle like-sign correlator (cos(0 Q + <fip — 20 c )). 
This indicates higher order terms plus the cross terms in same-sign correlation A{A 2 ) 
contribute very little and can be neglected. However, our opposite-sign correlation 
A(A + A_) differs from the opposite-sign correlator significantly. The A(A + A_) is 
significantly negative, while unlike-sign correlator is close to zero. This suggests that 
the higher order terms plus the cross terms are important. Note that the event-plane 
resolution correction for higher order terms may not be correct in this case. But 
the event-plane resolution effect can only smear out the correlations. Thus, the real 
magnitude of the A(A + A_) should be more negative with respect to the true reaction 
plane. See the discussion in section 13.5.11 

In order to make direct comparison to the published result, we analyze the charge 
asymmetry correlations within the entire TPC region. To avoid self-correlation, we 
randomly divide the event into two halves (sub-events) regardless of the 77 region. 
We use one sub-event to calculate the asymmetries and their correlations relative 
to the event-plane reconstructed from the other sub-event. We also calculate the 
three-particle correlators using particles a and from the sub-event for asymmetry 
calculation, and particle c from the sub-event for the event-plane reconstruction in 
the asymmetry calculation regardless of the charges. The results are shown in figure 
14.141 where the solid points are our LR — UD asymmetry correlations; the open points 
are our calculated the three-particle correlations; and the solid lines are the published 
data from [20]. Our calculated three-particle correlators match the published data 
very well, so we have confidence in the correctness of our analysis code. We also 
observe the qualitatively similar results to what we show in figure I4.13[ namely that 
the same-sign correlation is consistent with the like-sign three-particle correlator, and 
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Figure 4.14. This figure shows the same comparison as figure H~T3l but 
with the results obtained by randomly dividing the event into two equal 
sub groups of tracks within the full t] range of the TPC (—1 < rj < 1). 
The asymmetries and correlators are then calculated from one sub group 
with respect to the event-plane reconstructed from the other sub group. 
The solid curves are the published three-particle correlator values [20] 
superimposed upon our results. Error bars are statistical only. 



opposite-sign correlation is significantly more negative than the unlike-sign three- 
particle correlator. So we conclude that the high order terms and the cross terms in 
the UD — LR Fourier expansion are negligible for same-sign correlations, however, 
they are significant for opposite-sign correlations. 

4.5.2 Interpretation of Charge Asymmetry Correlations 

Our charge asymmetry correlations provide an opportunity to study the event- 
plane dependence of the charge separation effect. The published three-particle corre- 
lators assume that the higher order terms and cross terms are negligible. However, 
we showed in the previous section that the higher order terms and the cross terms 
contribute significantly to the opposite-sign correlations. Furthermore, our study can 
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provide more information on the charge asymmetry correlations in the in-plane and 
out-of-plane regions separately. In this section, we compare the CME/LPV expecta- 
tions to our results, and try to interpret what data tells us. 

First, we check the CME/LPV charge combination expectation. The CME/LPV 
would produce additional same-sign pairs in the up and down hemispheres. Such 
additional correlation will result in wider distributions (larger variances) of both the 
positive and negative charge multiplicities in the same hemisphere. Hence we expect 
larger dynamical correlation 5 (A 2 ) in UD than that in LR direction. On the other 
hand, the CME/LPV will produce back-to-back opposite-sign pairs in UD direction, 
which gives anti-correlation between the positive charge and negative charge parti- 
cle pairs in the UD direction. The result is a smaller covariance (A + AJ)ud than 
(A + A_) LR . 

Our data show that the same-sign dynamic correlation satisfy 5(A 2 JD ) > 5{A 2 LR ) 
of all centralities in Au+Au 200 GeV collisions (figure 14. 2p . This is qualitatively 
consistent with CME/LPV expectation because the multiplicity distribution is wider 
in UD direction than LR direction. However, the same-sign correlations are negative 
in mid-central to central collisions regardless of whether they are measured in-plane 
or out-of-plane. Negative variances suggest that the same-sign charged pairs are 
preferentially emitted back-to-back, in another words more symmetrically distributed 
between hemispheres. This is contrary to the expectation from the CME/LPV that 
same-sign charged pairs tend to be emitted in the same direction. 

The opposite-sign correlations show (A + A-)ud > (A + A—)lr in all centralities 
except for the most peripheral collisions, which suggests the opposite-sign particle 
pairs are likely emitted in the same hemisphere, and more strongly in UD than LR 
direction. This contradicts the naive expectation from CME alone that opposite-sign 
charged pairs are emitted back-to-back. The very peripheral collisions of Au+Au 
collision data are consistent with CME/LPV, however the same effect is also observed 
in d+Au collisions, where CME/LPV is not expected. The peripheral result can 
be explained by non-flow effect in low multiplicity events. When we divided an 
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event into two sub-events, the low multiplicity collisions are more sensitive to non- 
flow correlation, such as di-jet. The non-flow correlation will become important and 
dominate the reconstructed event-plane orientation with limited multiplicity. The 
event-plane is then preferentially sitting in the di-jet direction, and separating the 
event into two more equal halves in up and down hemispheres, which makes the di- 
jets lie in the LR hemispheres separately. This leads to a large fluctuation in the 
LR multiplicity, yielding large LR asymmetry. But for the UD hemispheres, they 
are divided more or less symmetrically by the event-plane axis, which has smaller 
asymmetries of the opposite-sign. Thus, the (A + A^) UD < (A + A—) LR for peripheral 
collisions is reasonable. 

Secondly, we look at the magnitude of the charge asymmetry correlation. There 
is little theoretical guidance to the quantitative magnitude of the charge separation 
itself and the charge correlations. However, there are model estimates suggesting a 
few percent charge asymmetry, which implies for our charge multiplicity asymmetry 
correlation magnitudes of the order of 10~ 4 to 10~ 3 . Considering the in-medium 
interaction and the multiplicity dilution effect, the estimated charge asymmetries 
could be reduced by one order in magnitude because the interaction with medium 
conserves parity and destroys only the correlations. This has been shown in |14j . 
Or, it could be even lower than 10 -6 as shown in [15J. Our results do indicate 
much larger correlations than these estimates, see figure 14.31 (divide the data points 
with Np ar t from 13. ip . Both the variance and covariance in mid-central collisions 
have the magnitudes around ~ 10~ 4 . There are also models arguing that the charge 
asymmetries and the magnitude can be explained by known QCD processes without 
invoking CME. Nevertheless, the estimated charge asymmetry correlations are, at 
least, a few orders of magnitude smaller than the measurements. Our measurements 
show strong correlations between same-sign and opposite-sign correlation in both U D 
and LR directions. 

Thirdly, we look at the pr dependence of the charge asymmetry correlations. Fig- 
ure |4]5] shows that the charge asymmetry UD — LR correlations grow with px for both 
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the same-sign and opposite-sign. While CME/LPV expects the charge separation to 
take place in low-pr region with p T < 1 GeV/c [TU], the data are inconsistent with 
that naive expectation. 

The wedge location dependence will be discussed in the following section 14.5.31 
Despite that both same-sign and opposite-sign correlations in figure H~3l are all pos- 
itive except the most peripheral collisions, the same-sign correlation A(A 2 ) is larger 
than the opposite-sign correlation A(A + A_). A yet unknown underlying background 
which is coupled with flow but exhibits no charge dependence could produce similar 
correlation. In that case, the correlation from underlying background could possibly 
lie in between our same-sign and opposite-sign correlation results. Thus, after sub- 
tracting the common background, the same-sign and opposite-sign correlations will 
have different sign, which is then consistent with CME/LPV expectation, assuming 
medium effect does not change correlation signs and the common background does 
not depend on charge combinations. Taking this idea, we define the charge separation 
observable (A) across the event-plane as 

A = A(A 2 ) - A(A + A^). (4.6) 

It is the correlation difference between the same-sign and opposite-sign UD — LR 
correlations. Equation 14.61 would quantify the charge separation effect if the reaction- 
plane dependent backgrounds were the same for the same-sign and opposite-sign 
correlations. A should be consistent with if there is no charge separation. A > 
would suggest a charge separation effect with same-sign charged particles emitted in 
the same direction (small angle correlation) and/or opposite-sign charged particles 
emitted in the opposite direction. A < would suggest the contrary situation with 
same-sign charge particles emitted in back-to-back direction and/or opposite-sign 
charged particles emitted in the same direction. We further test A in the following 
sections. 
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4.5.3 In-Plane or Out-of-Plane? 

If CME/LPV effect in the same-sign correlation could survive through the hydro- 
dynamic evolution to the final state, and its preferred direction is still along the sys- 
tem orbital angular momentum direction (magnetic field direction), then the charge 
separation A from equation 14.61 should depend on the wedge size and is possible to 



measure, as illustrated in figure 4.15(a) The smaller the wedge size, the stronger 



charge separation effect we should expect, for 

A(A0 w ) = A(Al^) - A(A + A.) A ^. (4.7) 

It is also possible that in-medium interaction and final state interaction will modify 
the preferred direction of the same-sign asymmetry correlation, so that A(A0 w ) may 
not necessarily increase with decreasing the wedge angle A0 W . For example, the 
elliptic flow pushes the correlated pairs into the in-plane direction, as illustrated in 



figure 4.15(b) 



The A(A0 w ) result is shown in figure KM of 20-40% Au+Au 200 GeV collisions. 
Charge separation decreases with decreasing the wedge size A0 W . The charge sepa- 
ration seems to disappear at very small wedge open angle in mid-central collisions. 
The same-sign and opposite-sign pair correlations have the maximum difference in the 
hemispheres. This may suggest that the effect of charge separation across the event- 
plane happens in the vicinity of the in-plane direction rather than the out-of-plane 



direction in mid-central collisions. Figure A. 19(d) shows the central 0-20% cen- 
trality charge separation wedge size dependence A(A0 w ), and figure A. 20(d) shows 
the peripheral 40-80% centrality A(A0 w ). Charge separation decreases with decreas- 
ing wedge size for all centralities. It even goes negative in central and peripheral 
collisions. 
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(a) Charge separation out-of-plane (b) Charge separation in-plane 



Figure 4.15. Schematic plot depicting the possible charge separation 
direction. The wedge size dependence of the difference between same-sign 
and opposite-sign UD — LR correlations (A(A0 w )) will increase with the 
decreasing of the wedge size, if the charge separation is out-of-plane as 
shown in panel (a). Otherwise, A(A0 w ) decreases with the decreasing of 
the wedge size, if the charge separation is in-plane as shown in panel (b). 
Our results favor in figure (b). 



4.5.4 Signal or Background? 

As discussed above, the charge separation along the system angular momentum 
direction (UD) is possibly due to CME, if the event-plane dependent common back- 
ground lies between the same-sign and opposite-sign correlation in figure 14.31 How- 
ever, we find inconsistencies with the naive CME/LPV expectation. Firstly, the 
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Figure 4.16. The wedge size dependence of the charge separation A(A0 w ) 
of Mid-central 20-40% Au+Au collisions. Error bars are statistical. 



observed same-sign correlation shows more back-to-back symmetry in both UD and 
LR direction; it is stronger in UD than LR. Secondly, the observed opposite-sign 
correlation shows that the opposite charged pairs are emitted in the same direction, 
and more strongly emitted in the UD than LR direction. Thirdly, the charge sep- 
aration seems to happen in the vicinity of the event-plane instead of out-of-plane. 
Fourthly, the asymmetry correlations increase with transverse momentum. Except 
the pt dependence, the results may be explained by CME/LPV if an unsubtracted 
background in these observables is charge independent and falls in between the same- 
and opposite-charge measurements. 

On the other hand, the assumption that the event-plane dependent background 
is charge independent may not be true. It is possible that the background for the 
same-sign and opposite-sign correlations is quite different due to different physics 
mechanisms. As argued in [40] by Pratt, charge conservation combined with collective 
flow naturally creates charge separation, with a difference between same-sign and 
opposite-sign correlations. It is already shown in figure I4T31 that the correlation UD — 
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LR differences as a function of centrality are qualitatively similar to the elliptic flow x>i 
centrality dependence [TS]. The correlations are peaked at medium central collisions, 
and drop close to zero in peripheral and most central collisions. Thus, we plot the 
charge separation A (defined in equation 14.61) against the average event-by-event 
anisotropy (v% bs ). The (v^ 38 ) value of each centrality is listed in table I3TT1 The result 
is shown in figure 14.171 The charge separation is scaled by number of participants 
Npart- We show the centrality bin numbers under the data points according to table 
13.11 The dashed line is a linear fit of the data points, which shows a very good 
agreement of the event shape dependence of the charge separation. The solid line is 
also a linear fit but with fixed intercept at zero. The result is qualitatively consistent 
with Pratt's suggestion that the charge separation is proportional to the elliptic flow 
in [3D], which is also suggested in [19] and consistent with the finding in STAR 
previous measurement in [TB]. The linear dependence of the charge separation effect 
on v^" 3 suggests that there may be a v^-dependent background in our measured 
charge separation. 

Recall the correlation dependence of vf" 3 in figure I4T71 The event-by-event anisotropy 
charge asymmetry correlations also suggest that the different behavior of the same- 
sign and opposite-sign correlations may be subject to different physics mechanism. 
The background may be due to a significant bulk correlation related to the event 
structure, which causes the background to be indeed different for same-sign and 
opposite-sign correlations. For proper interpretation of figure 14 .7\ it is important 
to understand the charge separation background and signal, and the cause of the 
possible signal. 

For the \ow-px v% dependent same-sign correlation, the increasing trend in UD 
direction suggests more and more back-to-back same-sign charged particle pairs are 
emitted in-plane while the event is more elongated in the in-plane direction, i.e. 
increasing v^ s . Equivalently, there might be fewer particles emitted out-of-plane 
with the increase of event-by-event anisotropy v^ s . It is also possible both are true. 
Thus the overall effect would be fewer back-to-back same-sign pairs emitted in out- 
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Figure 4.17. Charge separation A scaled by N part as a function of the 
measured average anisotropy (vf s ) in RUN IV 200 GeV Au+Au collisions. 
The centrality bin is labeled below each data point (see table I3TTT) . The 
particle pt range of 0.15 < px < 2.0 GeV/c is used for both asymmetry 
calculation and event-plane reconstruction. Error bars are statistical. 



of-plane direction than in-plane, causing the dynamic correlation in UD to increase 
with increasing v^ s . The decreasing trend in LR direction would be similar to that 
in UD direction. There are more particles emitted in the in-plane direction than 
out-of-plane direction with larger v^ s , thus more abundant back-to-back same-sign 
particle pairs in-plane than out-of-plane, which reduces the dynamical correlation in 
LR. 

For the opposite-sign correlations, the trend seems opposite to the same-sign cor- 
relations, but with much weaker effect. This may suggest different origins of the 
same-sign and opposite-sign particle pairs. The behaviors of S{A 2 ) and (A + A_) as 
a function of low-pr v^ s could also be a possible effect of cluster or resonance decay 
correlations overlaid with the elliptic flow [39] . 
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Figure 4.18. Panel (a): RUN IV Au+Au 200 GeV 20-40% centrality 
charge separation A as a function of low-pr event-by-event anisotropy 
v^ s . Panel (b): Event-plane resolution squared as a function of v^ s . The 
asymmetries and are calculated relative to the second order event- 
plane reconstructed from the other side of the TPC tracks. The particle 
Pt range of 0.15 < pt < 2.0 GeV/c is used for asymmetry, v^ s and 
event-plane reconstruction. Error bars are statistical. 
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We study the charge separation A as a function of the low-]>r event-by-event 
anisotropy vf" 3 to show the event shape dependence. In figure 4.18(a), we plot the 
charge separation between the same-sign and opposite-sign correlations A = A(A 2 ) — 
A(A + A_) as a function of the \ow-px event-by-event elliptic anisotropy v^ s . It is 
simply the difference between the same-sign and opposite-sign correlations shown in 



figure 4.7(c) Here v^ s is calculated on the event-by-event basis from particles used 
in asymmetry calculation relative to the event-plane reconstructed from the other 
half side of TPC tracks. Data are from RUN IV 20-40% centrality Au+Au 200 GeV 
collisions. For those events with large positive v^ s , the charge separation is large 
and positive, A(A 2 ) > A{A + A-), which is consistent with CME/LPV. However the 
charge separation flips sign with more negative v^ s , which appears inconsistent with 
CME/LPV expectation. 

One can argue that, at very negative v% bs , event-plane has very poor resolution, 
thus, the reconstructed EP is more orthogonal to the real reaction-plane rather than 
align with it. The UD and LR direction are then flipped if it is the real situation. So, 
the charge separation for those events is actually positive. The event-plane resolution 
does vary with the particle elliptic anisotropy vf" 3 . To show the dependence of event- 
plane resolution on v^ s , the square of the EP resolution from the half event is shown 



in figure 4.18(b) We obtain the resolution by randomly dividing the half event into 
two quarter-events noted by subscript a and b. Then we can reconstruct the event- 
plane from the quarter events to get ipEP,a and ipEP,b- Similarly, the half event-plane 
resolution can be assessed by e 2 EP = 2 (cos 2 (ipEP,a — ipEP,b))- Although the particles 
for event-plane reconstruction and event anisotropy calculation are from different 
phase space, we can still see the correlation between them in the figure. Note that at 
significant negative v^ s , the (cos2{ip ep,o, — ip EP,b)) turns negative, which suggests that 
the reconstructed event-plane does not reflect the true reaction-plane. It might be 
more likely orthogonal to the reaction-plane rather than aligned with it. This could 
mean that the U D and LR are indeed flipped. Therefore, the results at very negative 
v% bs are also consistent with CME. 
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As shown, the resolution squared at v^" 3 ~ is sizably positive, which means the 
event-plane reconstructed with an isotropic half event on the other side is considerably 
good. The charge separation A vanishes or slightly positive at v^ s ~ 0. For events 
with modest negative v^ s [y^" 1 > —0.1 for this particular centrality), the event-plane 
resolution from two sub-events method is well defined. Within —0.1 < v^ s < 0, 
the charge separation is negative with reasonable EP resolution. And the charge 
separation has a good linear dependence on the event-by-event anisotropy v % bs , which 
is not expected by CME/LPV. We know that if we integrate over all events within 
20-40% centrality bins, the average {v^" 3 ) is positive with the magnitude of a few 
percent, see table [37T1 The charge separation of the integrated result will be positive 
due to the linear dependence, but it doesn't necessarily mean the charge separation 
is caused by CME/LPV. The v^ s dependent charge separation result suggests that 
CME/LPV may not be the reason for the charge separation observed in the study. 

We show the charge separation vs vf" 1 results for central collision 0-20% centrality 
in figure IA.21| and for peripheral collisions 40-80% centrality in figure IA.22I The 
results are qualitatively similar to the mid-central results. 

It has been pointed out by several authors that there could be charge dependent 
physics background proportional to the event-by-event v^ s due to the net effect of 
particle intrinsic correlation and production anisotropy [16l[T8lll0]. The charge de- 
pendent background will cause the final state charge separation. We thus fit the 
charge separation A(v% bs ) to a linear polynomial in v^ s with the fitting range of 



—0.15 < vf" 3 < 0.15. The result is shown as red line in figure 4.18(a) We have 
A(>f s ) = (0.23 ±0.12) x 10" 4 + (3.0±0.2) x 10" 3 wf s . The slope could be a measure- 
ment of the particle intrinsic correlation strength. And the intercept could then be a 
more sensitive measurement of the charge separation. In other words, the measure- 
ment of CME/LPV has to take place in the phase space where the particle multiplic- 
ity distribution is more isotropic, so that the measured charge separation could be 
more sensitive to CME/LPV. To do so, we apply a cut requiring the event-by-event 
< 0.04 for the asymmetry calculation, and then plot the charge separation A 
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(b) Event-plane resolution of isotropic events 



Figure 4.19. Panel (a): The charge separation A scaled by N part as a 
function of centrality requiring the event isotropic condition \v% bs \ < 0.04. 
Panel (b): The event-plane resolution for those isotropic events. Data 
are from RUN IV 200 GeV Au+Au collisions, and particle pr range of 
0.15 < px < 2.0 GeV/c is used for asymmetry calculation, v^ s calculation 
and event-plane reconstruction. Error bars are statistical. 



as a function of centrality. The result is shown in figure 4.19(a), with the charge 
separation A scaled by N part . The charge separation is consistent with zero within 
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our present statistical precision, which suggests no substantial charge separation is 
observed in those approximately isotropic events of the measured particles. By taking 
the intercept value of the linear fit fuction of charge separation A = 0.23 x 10~ 4 at 
v^ s = 0., with twice of the uncertainty a A = 0.12 x 10~ 4 away from the intercept, we 
get the upper limit of the charge separation for the isotropical events in mid-central 
collisions, which is A = 4.7 x 10" 5 with 98% CL. 



The event-plane resolution for those nearly isotropic events is shown in figure 4.19(b) 



which shows that, with one spheric shape event on one side of the TPC, the event- 
plane reconstructed from the other side still has a good resolution. 

For consistency, we check the charge separation vs v^ s with one unit r\ gap result. 



It is shown in figure 4.20(a) and the corresponding event-plane resolution squared is 



shown in figure 4.20(b) With the large pseudo-rapidity gap, we further remove the 
short range bulk correlations from the soft particles. The charge separation fitted to a 
linear polynomial as a function of vf s gives (0.33 ±0.22) x 10~ 4 + (3.3 ±0.3) x 10" 3 ^ 6s 
with the fitting range of —0.15 < v^ s < 0.15. The linear dependence of the charge 
separation with event-by-event vf s is still seen after removing partially the short range 
correlations. Note for those isotropic events with v^ s ~ 0, the charge separation is 
close to zero with large errors due to the statistics. And resolution squared of the 
event-plane reconstructed from the half unit of pseudo-rapidity is reasonable positive 



at v% ~ as shown in figure 4.20(b) 



The central (0-20% centrality) and peripheral (40-80% centrality) collision results 
are shown in figure IA.23I and figure IA.24} which also show the charge separation is 
linearly dependent on the event-by-event v^ s . 

We can also study the charge separation using the first order event-plane recon- 
structed from the ZDC-SMD detector, which is not correlated to the TPC tracks. The 
charge separation with respect to the first order event-plane result is shown in figure 



4.21(a) along with the linear fit of the charge separation against v^" 1 also with respect 
to the first order event-plane within the range of —0.15 < v^ 11 < 0.15. We observe 



linear dependence of charge separation (0.03 ± 0.06) x 10 4 ± (3.6 ±0.1) x 10 3 t> 



3„,obs 
2 5 



115 



1.5 



x10" 



▲ A = A(A 2 )-A(A + A.) 

[0.33±0.22 + (33.1 2.9)x] 10 J 
X 2 /NDF=16.8/13 




-0.2 0.2 

<cos2((]>-\|f E p)> (p T <2GeV/c) 

(a) Charge separation vs v^ 3 with r\ gap 



CM 
</> 
O 

cS 



0.5- 



0- 



-0.5 


















••••• 










-0.2 C 


I 0.2 



(cos2((|)-^ E p)) (p T <2GeV/c) 

(b) EP resolution of 77 gap 



Figure 4.20. Panel (a): The charge separation A scaled by N part as a 
function of event-by-event anisotropy v^ s . The asymmetries and event- 
plane reconstruction are taken place in sub events with one unit pseudo- 
rapidity separation, —1.0 < 77 < —0.5 and 0.5 < 77 < 1.0. The charge 
separation is fitted to a linear polynomial as shown in red line. Panel (b): 
The event-plane resolution squared as a function of v^ s . Data are from 
RUN IV 200 GeV Au+Au collisions in 20-40% centrality, and the particle 
p T range of 0.15 < pr < 2.0 GeV/c is used for asymmetry calculation, 
v^ s calculation and event-plane reconstruction. Error bars are statistical. 
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Figure 4.21. Panel (a): The charge separation A scaled by N part as 
a function of event-by-event anisotropy v^ s . The asymmetries and v^ s 
are calculated from half TPC tracks of an event, with respect to the first 
order event-plane reconstructed from ZDC-SMD detectors. The charge 
separation is fitted to a linear polynomial as shown in red line. Panel (b): 
The first order event-plane resolution squared as a function of vf" 3 . Data 
are from RUN VII 200 GeV Au+Au collisions in 20-40% centrality, and 
the particle px range of 0.15 < p T < 2.0 GeV/c is used for asymmetry 
calculation and vf s calculation. Error bars are statistical. 
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which suggests the charge separation is indeed correlated with the event-by-event 
shape. 

At v^ s ~ 0, the charge separation is consistent with zero, indicating no charge 



separation effect for those spherical shaped events. Figure 4.21(b) shows the first 
order event-plane resolution squared as a function of v^ s . The large statistics of RUN 
VII data helps to limit the uncertainty while the first order event-plane resolution is 
relatively low compared to the second order event-plane. 

We also show the first order event-plane central (0-20% centrality) and peripheral 
(40-80% centrality) collision results in figure IA.25I and figure IA.26I They are quali- 
tatively consistent with the linear dependence and the intercepts are consistent with 
zero. 

Results for the top 2% most central collisions are shown in figure I4.22[ exhibiting 
charge separation that linearly depends on the event-by-event v^ s within the range 
of —0.1 < v% bs < 0.1, though the slope is smaller than the other three measurements 



above in figure 4.22(a) The event-plane resolution of the top 2% central collision 



is shown in 4.22(b) The linear dependence is qualitatively similar to the results 
presented earlier. 

To summarize, we have studied four cases of charge separation (A) as a function of 
event-by-event anisotropy {v^ s ). The results agree well with a linear dependence on 
the v^ s . The charge separation is larger with larger event elliptic shape v^ s , meaning 
the more elongated shape in in-plane direction, the more same-sign pairs going in 
the same direction and/or more opposite-sign pairs going back-to-back across the 
event-plane. The charge separation turns negative when v^ s goes negative, meaning 
the more elongated shape in out-of-plane direction, the more same-sign pairs going 
back-to-back and more opposite-sign pairs going in the same direction across the 
event-plane. Such linear effect is not expected from CME/LPV, thus it cannot be 
explained by CME/LPV alone. 
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Figure 4.22. Panel (a): The charge separation A scaled by N part as 
a function of event-by-event anisotropy t>2 6s of RUN IV ZDC triggered 
top 2% centrality 200 GeV Au+Au collisions. The asymmetries and v^ s 
are calculated from half TPC tracks of an event, with respect to the 
event-plane reconstructed from the other half TPC tracks. The charge 
separation is fitted to a linear polynomial as shown in red line. Panel (b): 
The top 2% centrality event event-plane resolution squared as a function 
of v^ s . The particle px range of 0.15 < pr < 2.0 GeV/c is used for 
asymmetry calculation and v^ s calculation. Error bars are statistical. 
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5. SUMMARY 

Motivated by the Chiral Magnetic Effect (CME), a possible signature of local parity 
violation (LPV) in heavy ion collisions, we have measured the charge multiplicity 
asymmetries and their correlations for both positively and negatively charged parti- 
cles in the directions across the event-plane (UD) and the plane perpendicular to the 
event-plane (LR). Each event is divided into two equally halves (sub-events) accord- 
ing to the pseudo-rapidity of the tracks within TPC (— 1 < rj < and < rj < 1). 
The event-plane is reconstructed from one sub-event, and the asymmetries are cal- 
culated in the other sub-event with respect to the event-plane in order to reduce 
self-correlation. The dynamical variances (5(A±)) and covariances (8(A + A_)) of pos- 
itive and negative charges are presented in UD and LR directions with Au+Au and 
d+Au at 200 GeV center of mass energy collisions. We use asymmetry correlation 
from LR direction as a null- reference to study the CME/LPV which is an effect 
supposedly along the system angular momentum direction, i.e. UD direction. 

As shown in figure H~2"| the same-sign dynamical variances S(A 2 ) are positive in 
d+Au collisions and peripheral Au+Au collisions both out-of-plane (UD) and in-plane 
(LR), and turn negative in mid-central and central collisions. The positive variances 
in d+Au and peripheral Au+Au collisions indicate a broadening in the multiplic- 
ity asymmetry distributions, which suggests the same-sign charged particle pairs are 
positively correlated, or in other words, preferentially emitted in the same direction 
(small angle pairs). On the other hand, in mid-central and central collisions, the 
variances are turning negative, which indicates narrowing in the multiplicity asym- 
metry distributions, suggesting the same-sign charged pairs are negatively correlated, 
and preferentially emitted in the back-to-back direction. The dynamical asymmetry 
variance 5(A 2 JD ) is larger than S(A 2 LR ) except the most peripheral bin. This indi- 
cates that there are more same-sign small angle pairs emitted in out-of-plane than 
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in-plane direction, equivalently, more same-sign back-to-back pairs in in-plane than 
out-of-plane. 

The opposite-sign covariances (A + A_) are largely positive in d+Au and all Au+Au 
collisions both in-plane and out-of-plane. It is an indication of strong positive correla- 
tion of opposite-sign charged particle pairs, which suggests the positively charged and 
negatively charged particles are preferentially emitted in the same direction regard- 
less of the reaction plane direction. There are more significant positive correlation 
in out-of-plane than in-plane direction, (A + AJ)ud > (A + A_) lr, for most centrali- 
zes except the most peripheral two bins, which is consistent with d+Au collisions. 
The additional positive asymmetry correlation indicates opposite-sign charged pairs 
in out-of-plane are more preferentially emitted in small angle than that in in-plane 
direction. 

By taking the difference between the correlations out-of-plane (UD) and in-plane 
(LR), we obtain asymmetry correlations UD — LR by subtracting the charge depen- 
dent correlation background which is not related to the reaction-plate as shown in 
figure POl The variance difference A{A 2 ) = 6(A 2 ID ) — 5(A 2 )lr and covariance differ- 
ence A(A + A^) = (A + A-)ud — (A+A-) lr show similar centrality dependence, and 
are all positive except most peripheral covariance bins. The UD — LR result shows 
inconsistence with CME/LPV for the opposite-sign correlation, which is naively ex- 
pected to be negative due to CME/LPV charge separation. The same-sign variance 
is consistent with the CME/LPV expectation that there is additional broadening of 
charge multiplicity asymmetry in out-of-plane direction compared to the in-plane di- 
rection. However, we find (section H~Tj) that the variances 5(A + A^) UD and 5(A + A_) L r. 
are negative in mid-central to central collisions, suggesting additional back-to-back 
same-sign pairs both in-plane and out-of-plane, which is inconsistent with CME/LPV 
expectation. 

The transverse momentum dependence figure WM shows that, the charge asym- 
metry dynamical variances 5{A 2 )ud and 5{A 2 )lr decrease with px- They are all 
positive at low p?, then decrease below zero for dynamical variances and covariances 
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above 1 GeV/c. Data indicates the small angle same-sign correlations are stronger in 
low p T , suggesting a bulk effect. The charge asymmetry covariances (A + A_) UD and 
(A + A^) LR are all positively correlated and increase with p T , and (A + A_) UD grows 
faster than (A + A^) LR with p T . This suggests, for the entire px range, opposite-sign 
particle pairs are always positively correlated, i.e. positive and negative charges are 
preferentially emitted in same direction, and with stronger effect in higher p?. 

STAR has previous measurements of three-particle correlators which were consis- 
tent with CME/LPV expectation together with medium interaction. However, the 
possible physics backgrounds which are related to reaction-plane have not been ad- 
dressed thoroughly. We investigated the connection between our charge asymmetry 
correlation measurements and the three-particle correlator measurements, and showed 
that the differences are due to higher orders and cross terms in the Fourier expansion 
of asymmetry correlations as a function of azimuth, especially in the opposite-sign 
covariances, shown in figure 14.131 and figure 14. 141 

To investigate charge separation, we study the wedge size and wedge location 
dependence of the charge asymmetry correlations. The charge asymmetry correlations 
as a function of the wedge size is shown in figure 14. 1 11 Both the dynamical variances 
and covariances increase with decreasing wedge size, which suggests the reaction- 
plane dependent charge asymmetry correlations are more likely local in azimuth. The 
UD — LR measurements show the difference both A(A + A_) and A(A 2 ) increase with 
decrease wedge size. However, the difference between same-sign and opposite-sign 
correlations vanish for small wedge size. 

A common unknown background for the same-sign and opposite-sign correlations 
could possibly lie in between the correlations. Then the difference between variance 
and covariance A = A(A 2 ) — A(A + A^) , might be sensitive to the charge separation 
induced by CME/LPV. We find that the charge separation A^ w decreases with 
decreasing of wedge size, figure 14. 16} which suggests the charge separation across the 
event-plane happens in the vicinity of the in-plane rather than out-of-plane direction. 



122 



We have also studied the charge asymmetry correlations as a function of the event- 
by-event anisotropy v^ s of the measured asymmetry particles in figure 14.61 and figure 
14.71 The UD — LR correlations show little dependence on high-p^ v^ s . However, the 
UD — LR of the variances A(A 2 ) increases strongly with low-pr v 2 bs > while the covari- 
ance A(A + A_) decreases slowly with low-p? v^ s . Also the variance and covariance 
intercept at the same positive value at v^ s « 0. 

We reported the charge separation A as a function of v^ s with different event- 
plane reconstruction methods and event selection cuts. The results show a robust 
linear dependence on v^ s . For all the instances, the charge separation across the 
event-plane A is qualitatively consistent with zero for those events with least event- 
by-event anisotropy {v^ ~ 0). For the most isotropic events (|t>2 bs | < 0.04) where 
such backgrounds may vanish, we find the charge separation effect is consistent with 
zero as shown in figure 14.191 

We show the charge separation A as a function of the average vf" 1 , (v^ s ), in figure 
14.171 A very good linear dependence is observed over the event shape, which suggests 
an event anisotropy dependent charge separation. 

It is also possible the physics backgrounds are different for same-sign and opposite- 
sign UD — LR correlations. The linear dependence of charge separation on event-by- 
event anisotropy v^ s suggests that the intrinsic particle correlation and event shape 
can play an important role in the charge asymmetry correlations. Also we suggest that 
alternative contributions can naturally create such effect as well. By taking the events 
with minimum anisotropy (v^ ~ 0), we found the charge separation magnitude to 
be A = (2.3 ± 1.2) x 10~ 5 for mid-central collisions (20-40% centrality), from figure 
14.181 Then, we conclude that CME/LPV can't explain the charge separation alone, 
and we provide the upper limit with current measurements with charge separation 
A = 4.7 x 10" 5 with 98% CL. 
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Figure A.l. Single track azimuthal angle distributions for all centralities before any acceptance corrections. Red 
and blue lines are for positive charged particles from 77 > and 77 < regions, and green and black lines are 
negative charged particles from 77 > and 77 < regions respectively. 
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Figure A. 2. Single track azimuthal angle distributions for all centralities after the acceptance corrections. Red 
and blue lines are for positive charged particles from 77 > and 77 < regions, and green and black lines are 
negative charged particles from 77 > and 77 < regions respectively. 
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Figure A. 3. Reconstructed second order event-plane azimuthal distributions for RUN IV Au + Au 200 GeV 
collisions in all centralities. The event-plane is reconstructed from charged particles within 0.15 < < 2 GeV ' jc 
from 7] < (red) and rj > (blue) separately. Error bars are statistical only. 
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Figure A. 4. Asymmetry correlations: panel (a) (A+ ud ) v >q, panel 
(b) (A 2 +MD ) V<0 , panel (c) (A+AJ) UD , panel (d) (A + A_) UD , panel (e) 
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from one side of the TPC, and the asymmetry correlations are calculate in 
the same pt range but from the other side of the TPC. 
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Figure A. 7. Charge multiplicity asymmetry correlations (A 2 ) as a function of the event-plane resolution €ep in 
all centralities. The solid lines are linear fits to the data. Error bars are statistical. 
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Figure A. 8. Charge multiplicity asymmetry correlations (A + AJ) as a function of the event-plane resolution €ep 
in all centralities. The solid lines are linear fits to the data. Error bars are statistical. 
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Figure A. 9. Charge multiplicity asymmetry correlations differences between UD and LR as a function of the 
event-plane resolution e^p in all centralities. The solid lines are linear fits to the data. The dashed lines are 
linear fits with fixed zero intercept at €ep = 0. Error bars are statistical. 
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Figure A. 10. Reconstructed first order event-plane azimuthal distributions for RUN VII Au+Au 200 GeV 
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Figure A. 11. Panel (a): Statistical fluctuation and detector effects in 
charge asymmetry variances scaled by the number of participants N part from 
east-side of the TPC, 77 < region, with respect to the EP reconstructed 
from west-side of the TPC, 77 > region. The black curve shows the 
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fluctuation plus detector effects (A 2 + LR stat+det ) are shown in green crosses 
with scramble method, and red circles with flipping- ir method. Panel (b): 
Same as (a) but for 77 > region with EP reconstructed from rj < region. 
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Figure A. 12. Panel (a): Statistical fluctuation and detector effects in 
charge asymmetry variances scaled by the number of participants N part from 
east-side of the TPC, 77 < region, with respect to the EP reconstructed 
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Figure A. 13. Asymmetry correlations vs event-by-event of RUN IV 
200 GeV Au+Au central 0-20% collisions. Error bars are statistical only. 
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Figure A. 14. Asymmetry correlations vs event-by-event of RUN IV 
200 GeV Au+Au peripheral 40-80% collisions. Error bars are statistical 
only. 
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Figure A. 15. Asymmetry correlations vs event-by-event v^ s of RUN IV 
200 GeV Au+Au central 0-20% collisions. The particles used for asymmetry 
calculation and event-plane reconstruction are divided by pseudo-rapidity 
— 1.0 < t] < —0.5 and 0.5 < t] < 1.0. Error bars are statistical only. 
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Figure A. 16. Asymmetry correlations vs event-by-event v^ s of RUN IV 
200 GeV Au+Au peripheral 40-80% collisions. The particles used for asym- 
metry calculation and event-plane reconstruction are divided by pseudo- 
rapidity — 1.0 < 7] < —0.5 and 0.5 < rj < 1.0. Error bars are statistical 
only. 
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Figure A. 17. Asymmetry correlations vs event-by-event v^ s of RUN VII 
200 GeV Au+Au central 0-20% collisions. The particles used for asymmetry 
calculation are from half side of the TPC with respect to the first order 
event-plane reconstructed from ZDC-SMD. Error bars are statistical only. 
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Figure A. 18. Asymmetry correlations vs event-by-event v^ s of RUN VII 
200 GeV Au+Au peripheral 40-80% collisions. The particles used for asym- 
metry calculation are from half side of the TPC with respect to the first 
order event-plane reconstructed from ZDC-SMD. Error bars are statistical 
only. 



20 40 60 80 
Wedge size (degree) 



qI 

100 T3 20 40 60 80 

Wedge size (degree) 



100 



(a) Asymmetry correlations vs wedge size (b) UD — LR correlations vs wedge size 



o 
o 

E 

>> 

< 



,x10" ; 



0- 



, , , , 


I ' ' 


I l l 










• 


(AAV* J 




■ 










■ 

. I . . 


I . . 


i ... i ... i ... 



20 40 60 80 
Wedge location ^ (degree) 



100 




20 40 60 80 
Wedge size A^ w (degree) 



100 



(c) Asymmetry correlations vs wedge location 



(d) Charge separation A(A0 w ) 



Figure A. 19. The wedge size dependence of asymmetry correlations in 
panel (a) and their differences between UD and LR correlations A(A 2 A ^) 
and A(A|_A_)a</> w in panel (b). Wedge location dependence is shown in 
panel (c) with opening angle of 30° (A0 W = 15°). Charge separation 
A(A0 w ) as a function of the wedge size panel (d). Data are from 0-20% 
centrality RUN IV 200 GeV Au+Au collisions. Error bars are statistical. 
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Figure A. 20. The wedge size dependence of asymmetry correlations in 
panel (a) and their differences between UD and LR correlations A(A 2 A ^) 
and A(A|_A_)a</> w in panel (b). Wedge location dependence is shown in 
panel (c) with opening angle of 30° (A0 W = 15°). Charge separation 
A(A0 w ) as a function of the wedge size panel (d). Data are from 40-80% 
centrality RUN IV 200 GeV Au+Au collisions. Error bars are statistical. 
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Figure A.21. Panel (a): RUN IV Au+Au 200 GeV 0-20% centrality 
charge separation A as a function of low-pr event-by-event anisotropy v^ s . 
Panel (b): Event-plane resolution squared e 2 as a function of vf"*. The 
asymmetries and vf" 3 are calculated relative to the second order event- 
plane reconstructed from the other side of the TPC tracks. The particle px 
range of 0.15 < px < 2.0 GeV/c is used for asymmetry, v^ s and event-plane 
reconstruction. Error bars are statistical. 
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Figure A.22. Panel (a): RUN IV Au+Au 200 GeV 40-80% centrality 
charge separation A as a function of \ow~pt event-by-event anisotropy v^ s . 
Panel (b): Event-plane resolution squared e 2 as a function of vf"*. The 
asymmetries and vf" 3 are calculated relative to the second order event- 
plane reconstructed from the other side of the TPC tracks. The particle px 
range of 0.15 < pr < 2.0 GeV/c is used for asymmetry, v^ s and event-plane 
reconstruction. Error bars are statistical. 
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Figure A. 23. Panel (a): The charge separation A scaled by N part as a 
function of event-by-event anisotropy v^ s . The asymmetries and event- 
plane reconstruction are taken place in sub events with one unit pseudo- 
rapidity separation, —1.0 < rj < —0.5 and 0.5 < rj < 1.0. The charge 
separation is fitted to a linear polynomial as shown in red line. Panel (b): 
The event-plane resolution squared as a function of v^ s . Data are from 
RUN IV 200 GeV Au+Au collisions in 0-20% centrality, and the particle 
Pt range of 0.15 < px < 2.0 GeV/c is used for asymmetry calculation, 
calculation and event-plane reconstruction. Error bars are statistical. 
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Figure A. 24. Panel (a): The charge separation A scaled by N part as a 
function of event-by-event anisotropy v^ s . The asymmetries and event- 
plane reconstruction are taken place in sub events with one unit pseudo- 
rapidity separation, —1.0 < rj < —0.5 and 0.5 < rj < 1.0. The charge 
separation is fitted to a linear polynomial as shown in red line. Panel (b): 
The event-plane resolution squared as a function of v^ s . Data are from 
RUN IV 200 GeV Au+Au collisions in 40-80% centrality, and the particle 
Pt range of 0.15 < pt < 2.0 GeV/c is used for asymmetry calculation, 
calculation and event-plane reconstruction. Error bars are statistical. 
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Figure A. 25. Panel (a): The charge separation A scaled by N part as 
a function of event-by-event anisotropy v^ s . The asymmetries and u£ 6s 
are calculated from half TPC tracks of an event, with respect to the first 
order event-plane reconstructed from ZDC-SMD detectors. The charge 
separation is fitted to a linear polynomial as shown in red line. Panel (b): 
The first order event-plane resolution squared as a function of v^ s . Data 
are from RUN VII 200 GeV Au+Au collisions in 0-20% centrality, and 
the particle px range of 0.15 < pt < 2.0 GeV/c is used for asymmetry 
calculation and v^ s calculation. Error bars are statistical. 
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Figure A. 26. Panel (a): The charge separation A scaled by N part as 
a function of event-by-event anisotropy v^ s . The asymmetries and v^ s 
are calculated from half TPC tracks of an event, with respect to the first 
order event-plane reconstructed from ZDC-SMD detectors. The charge 
separation is fitted to a linear polynomial as shown in red line. Panel (b): 
The first order event-plane resolution squared as a function of uJj ■ Data 
are from RUN VII 200 GeV Au+Au collisions in 40-80% centrality, and 
the particle px range of 0.15 < pt < 2.0 GeV/c is used for asymmetry 
calculation and v^ s calculation. Error bars are statistical. 
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